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Preface 


''A 

The  unique  properties  of  CaAs  make  it  possible  to  construct  inte¬ 
grated  circuit  devices  that  are  impossible  in  Si.  The  Air  Force  Avionics 
Laboratory/AADR  has  been  developing  this  technology  for  a  number  of 
years.  The  difficulty  of  introducing  dopants  by  diffusion  has  lead  ion 
implantation  to  play  an  increasing  role  in  the  fabrication  process.  The 
present  production  technique  for  high  performance  devices  is  to  fabricate 
large  quantities  and  select  those  few  that  meet  the  desired  specifica¬ 
tions.  Having  a  nondestructive  technique  that  can  be  used  to  charact¬ 
erize  the  implantation  process  during  fabrication  of  the  device  so  as  to 
reject  faulty  device  structures  can  save  valuable  time  as  well  as  money. 
Depth-resolved  cathodoluminescence  is  a  process  that  can  be  used  for  this 
purpose.  This  research  develops  and  verifies  a  model  of  cathodolumines¬ 
cence  in  ion  implanted  GaAs.  This  model  can  now  be  used  as  a  tool  for 
further  study  of  ion  implanted  GaAs.  This  is  the  first  step  in  devel¬ 
oping  cathodoluminescence  as  a  tool  for  deducing  the  shape  of  the  ion 
Implanted  depth  profile  in  semiconductor  materials. 

I  would  like  to  thank  my  advisor.  Dr.  R.  L.  Hengehold,  for  sug¬ 
gesting  the  investigation  of  this  area.  This  advice,  support,  and 
prodding  throughout  this  effort  are  the  primary  reasons  this  dissertation 
was  ever  completed.  1  owe  a  special  debt  to  Drs.  Y.  S.  Park,  B.  J. 
Pierce,  T.  Luke,  P.  E.  Nielsen,  and  J.  Jones,  Jr.  who  came  to  my  rescue 
at  several  critical  times  during  the  analysis.  I  am  also  grateful  for 
the  timely  and  skilled  assistance  given  me  by  Jim  Miskimen,  Ron  Gabriel 
and  George  Gergal  of  the  AFIT  physics  laboratory  staff  and  for  the  superb 
typing  by  Jill  Rueger. 
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I  would  also  like  to  give  a  very  special  thanks  to  my  wife,  Char¬ 
lene,  and  children.  Derrick  and  Catherine.  They  sacrificed  much  more 
than  I  in  order  to  obtain  this  goal. 
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Abstract 

Depth-resolved  cathodoluminescence  was  used  to  study  Ion  Implanted 
GaAs.  This  was  done  in  a  two  step  process.  First  a  model  for  the 
luminescence  from  ion  implanted  GaAs  was  developed.  This  model  Includes 
a  detailed  Monte  Carlo  simulation  of  electrons  penetrating  into  GaAs. 

The  result  of  this  calculation  is  a  prediction  of  luminescent  intensity 
as  a  function  of  electron  beam  energy.  This  calculation  differs  from 
preceeding  ones  in  that  the  ion  implanted  profile  is  specifically  in¬ 
cluded  in  the  analysis.  Second,  the  model  was  validated  by  a  suitable 
experiment  that  confirms  the  theoretical  predictions. 

The  luminescence  measurements  were  made  on  Mg  ions  implanted  into 
epitaxial  GaAs.  Mg  was  chosen  since  it  had  not  been  previously  studied 
in  epitaxial  GaAs. 

The  results  of  this  research  include  the  identification  of  many  of 
the  lines  in  the  Mg  implanted  GaAs  spectrum  from  1.32  ev  to  1.52  ev,  a 
study  of  the  effect  of  changing  the  current  of  the  electron  beam  on  the 
spectra  and  a  comparison  of  the  theoretical  and  experimental  luminescence 
curves.  These  results  show  that  changing  the  current  density  at  the 
sample  surface  dramatically  changes  the  spectra  and  that  good  agreement 
exists  between  the  theoretical  and  experimental  luminescence  curves. 
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1.  Introduction 


Overview 

There  is  an  increasing  emphasis  on  the  use  of  GaAs  as  a  follow-on  to 
silicon  in  electronic  devices.  This  lias  occurred  for  two  reasons.  First, 
some  of  the  properties  of  CuAs  such  as  a  higher  mobility  and  larger  band 
gap  allow  important  improvements  (higher  frequency  and  temperature  opera¬ 
tion)  in  the  performance  of  classical  devices.  Second,  some  of  the  char¬ 
acteristic  features  of  CaAs  give  rise  to  physical  phenomena  such  as  a  high 
frequency  instability  (Gunn  effect)  and  light  emission  in  the  visible  and 
IR  by  carrier  recombination.  These  phenomena  have  recently  been  used  in 
such  types  of  devices  as  microwave  oscillators  and  amplifiers,  light  emit¬ 
ting  diodes,  and  lasers.  The  difficulty  of  introducing  dopants  into  GaAs 
by  diffusion  and  a  desire  to  use  shallower  device  structures  has  led  to  an 
increasing  role  for  ion  implantation  as  the  means  for  introducing  these 
dopants  into  GaAs. 

This  work  is  important  to  the  U.S.  Air  Force  since  it  is  involved  in 
the  use  of  GaAs  as  a  follow-on  to  silicon  for  semiconductor  components.  The 
primary  interest  is  in  using  GaAs  in  microwave  and  optoelectronic  devices. 
Having  a  nondestructive  technique  that  can  be  used  to  characterize  the 
semiconductor  device  "in  situ"  before  further  processing  can  save  the  Air 
Force  valuable  time  as  well  as  money. 

in  order  to  improve  the  understanding  of  the  basic  implantation  pro¬ 
cess  and  to  predict  the  effect  on  device  behavior  of  variations  in  the 
implantation  process,  several  diagnostic  techniques  have  been  developed.  In 
general,  except  for  electrical  measurements  using  the  Hall  technique. 
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radio-tracer  sectioning,  and  static  and  transient  C-V  measurements,  these 
techniques  use  either  photons,  ions,  or  electrons  as  a  probe  and  lumines¬ 
cence,  x-rays,  electrons,  gamma  rays,  or  ions  as  the  profile  sensing 
mechanism.  They  include  Rutherford  scattering,  ESCA  (Electron  Spectro¬ 
scopy  for  Chemical  Analysis),  Auger  spectroscopy,  SIMS  (Secondary  Ion 
Mass  Spectroscopy),  and  PRP  (Proton  Resonance  Profiling)  just  to  name  a 
few.  Reference  (11)  explains  how  some  of  these  techniques  can  be  used. 
All  of  these  techniques  suffer  from  one  or  more  of  the  following  limita¬ 
tions  : 

1.  They  destroy  the  sample  either  by  bombardment  with  high  energy 
particles  or  successive  layer  removal  in  order  to  obtain  depth  informa¬ 
tion  about  the  impurity  implant. 

2.  The  technique  requires  a  p-n  junction  be  formed  and  several 
samples  be  tested  for  one  profile. 

3.  Suitable  reactions  are  necessary  for  the  implanted  species. 

4.  The  profiling  process  requires  one  or  more  days  to  complete. 

5.  The  equipment  for  the  analysis  is  very  expensive. 

A  technique  that  can  be  used  for  profiling  which  eliminates  the 
above  problems  would  be  extremely  useful.  Such  a  technique  could  allow 
routine  testing  of  device  structures  in  an  early  fabrication  phase,  thus 
avoiding  further  processing  on  defective  ones.  This  would  save  consid¬ 
erable  time  and  money  in  fabricating  circuits  from  GaAs.  One  such  pro¬ 
cess  is  depth-resolved  cathodoluminescence.  This  process  uses  a  beam  of 
electrons  to  excite  carriers  in  the  target  that  then  recombine  radia- 
tlvely.  fly  using  low  current  densities  and  low  beam  energies  (20  keV  or 
less)  the  sample  is  not  altered.  The  test  equipment  necessary  for  the 
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technique  is  not  very  expensive  compared  to  some  of  the  other  profiling 
techniques  listed  above.  Depth-resolved  cathodoluminescence  does  not 
require  a  p-n  junction  be  formed  or  more  than  one  sample  tested  in  order 
to  obtain  the  implant  profile  versus  depth  into  the  sample.  While  there 
is  a  requirement  for  the  implant  species  to  be  optically  active  in  order 
to  be  used  for  cathodoluminescence,  this  is  not  a  practical  limitation 
since  the  common  implant  species  in  GaAs  are  optically  active  after 
annealing.  (Only  if  the  profile  is  desired  before  optical  activation  of 
the  sample  (annealing)  takes  place  is  cathodoluminescence  limited.) 
Lastly,  the  profiles  can  be  done  extremely  fast  with  the  proper  test  set¬ 
up.  By  automatically  sweeping  the  electron  beam  voltage  and  processing 
the  resulting  data  on  an  online  computer,  the  profile  can  be  obtained  in 
minutes. 

Cathodoluminescence  on  GaAs  has  been  performed  by  many  authors,  some 

of  which  are  listed  in  references  (1-6,  11-16).  Especially  Dumoulin^ 

observed  changes  that  indicated  that  cathodoluminescence  could  be  used  as 

a  technique  to  obtain  depth  profiles  of  the  dopant  in  ion  implanted  GaAs. 

His  results  were  qualitative,  however  quantitative  results  have  also  been 
(7  9  is  19  21) 

attempted  ’  ’  ’  ’  .  None  of  these  quantitative  efforts  have 

been  concerned  with  the  cathodoluminescence  from  ion  implanted  semicon¬ 
ductors.  This  current  study  tries  to  establish  a  quantitative  basis  for 
cathodoluminescence  from  such  implanted  semiconductors.  This  is  the 
first  step  in  developing  cathodoluminescence  as  a  tool  for  deducing  the 
shape  of  the  ion  implanted  depth  profile  in  semiconductor  materials. 
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Summary  of  the  Research 

The  objective  of  this  research  is  to  develop  and  verity  a  model 
describing  cathodoluminescence  in  ion  implanted  GaAs.  This  was  accom¬ 
plished  in  a  two  step  process-  First,  a  model  for  the  luminescence  from 
ion  implanted  GaAs  was  developed.  This  model  differs  from  preceding  ones 
in  that  the  implanted  ion  distribution  is  specifically  included  in  the 
analysis.  The  result  of  this  portion  of  the  research  was  used  to  predict 
the  luminescent  intensity  as  a  function  of  the  electron  beam  energy  for  a 
specific  ion  implanted  profile.  Second,  the  model  was  validated  by  a 
suitable  experiment  that  confirms  the  theoretical  predictions.  This 
constitutes  the  forward  problem,  that  is,  predicting  the  luminescence 
assuming  a  known  implant  profile.  The  inverse  problem  of  predicting  the 
implant  profile  versus  depth  given  the  luminescence  curve  has  not  been 
attempted  and  constitutes  an  area  for  future  work. 

The  luminescence  measurements  were  made  on  Mg  ions  implanted  in  epi¬ 
taxial  GaAs.  Epitaxial  GaAs  was  used  since  it  is  the  highest  purity  GaAs 
available  and  hence  there  are  fewer  transitions  competing  with  the  Mg 
line.  Mg  was  chosen  since  it  has  not  been  previously  studied  in  GaAs 
epilayers  and  since  theoretical  studies of  Mg  ion  penetration  into 
GaAs  indicate  it  should  penetrate  to  a  depth  comparable  to  that  of  the 
electron  beams  being  used. 

The  results  of  this  research  include  the  identification  of  many  of 


the  lines  in  the  Mg  implanted  GaAs  spectrum  from  1.32  eV  to  1.52  eV,  a 
study  of  the  effect  of  changing  the  current  density  of  the  electron  beam 
on  the  spectra  and  a  comparison  of  the  theoretical  and  experimental 
luminescence  curves.  These  results  show  that  changing  the  current  density 


at  the  sample  surface  Jr  imatically  changes  the  spectra,  and  that  good 
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agreement  exists  between  the  theoretical  and  experimental  luminescence 
curves . 

The  principal  conclusion  of  this  research  is  that  the  model  presented 
he*"e  of  the  cathodol uml nescence  experiment  can  be  used  to  accurately  pre¬ 
dict  the  relative  intensity  of  the  ion  implanted  line  in  epitaxial  CaAs. 

Organization  of  Dissertation 

The  remaining  chapters  of  Litis  dissertation  describe  the  model,  the 
experiment,  and  present  results,  conclusions,  and  recommendations.  In 
Chapter  II  the  analytical  model  for  cathodol uminescence  is  developed. 

This  model  consists  of  three  parts.  The  first  part  is  a  Monte  Carlo 
simulation  of  electrons  penetrating  into  CaAs.  This  simulation  provides 
the  information  necessary  to  determine  the  generation  function  for  elec¬ 
tron-hole  pairs.  The  electron-hole  pairs  diffuse  through  the  GaAs  before 
recombining.  Thus  the  development  and  solution  of  the  appropriate  dif¬ 
fusion  equation  is  covered  in  the  second  section  of  Chapter  II.  The 
third  section  contains  a  discussion  ot  the  equation  used  to  calculate  the 
luminescence  versus  beam  voltage,  L(V),  curves.  This  equation  considers 
the  shape  of  the  impLanted  ion  profile,  the  excess  carrier  density  pre¬ 
dicted  by  the  solution  to  the  diffusion  equation,  and  the  absorption 
coefficient  of  the  emitted  radiation.  In  addition,  all  L(V)  curves 
derived  from  the  ion  Implant  are  normalized  to  the  L(V)  curve  for  a 
uniform  Implant.  This  normalization  removes  variations  in  the  experi¬ 
mental  data  that  might  be  expected  to  occur  from  run  to  run.  This  makes 
the  comparison  of  experimental  and  analytical  results  easier.  Next,  a 


Study  of  how  the  L(V)  curves  vary  as  the  diffusion  length,  surface 
recombination  velocity,  absorption  coefficient  and  impurity  profile  are 
changed.  The  chapter  concludes  witli  a  summary  of  how  the  calculation  is 
performed. 

In  Chapter  ill  the  experimental  system  and  the  procedures  used  to 
verify  the  model  are  described.  The  experimental  system  necessary  to 
verify  the  model  requires  certain  features  in  order  to  get  any  kind  of 
consistent  results.  These  features  include:  an  electron  beam  which  can 
be  easily  and  repeatedly  controlled,  a  beam  current  density  which  must  be 
uniform  across  tike  sample,  and  an  optical  system  which  should  be  insensi¬ 
tive  to  small  variations  in  position.  In  the  first  section  a  description 
is  given  ut  an  existing  cathodolumineseence  system  and  how  it  was  modi¬ 
fied  to  meet  these  particular  requirements.  The  following  section  contains 
a  review  of  the  previous  work  with  Mg  doped  and  implanted  GaAs.  All  pre¬ 
vious  work  has  been  with  substrate  quality  GaAs .  This  is  the  first 
report  of  extensive  measurements  on  Mg  implanted  epitaxial  GaAs .  The 
linal  section  of  Chapter  III  contains  a  list  of  the  samples  available  for 
testing,  how  they  were  made  and  how  they  were  processed  before  testing. 

The  resuLts  obtained  using  the  system  described  in  Chapter  111  and 
the  comparison  of  analytical  and  experimental  L(V)  curves  comprise  Chapter 
IV.  First,  the  spectrum  from  Mg  Ion  Implanted  GaAs  is  analyzed.  Various 
lines  between  1.3  and  1.5  eV  are  Identified.  It  is  shown  how  depth 
resolved  cathodolumineseence  can  be  used  to  identify  lines  by  determining 
where  the  concentration  of  the  Impurity  is  maximum.  Next,  a  study  is 
conducted  of  the  spectral  variation  with  current  density.  A  possible 
explanation  for  this  variation  is  given.  The  reason  that  this  variation 
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Is  Important  to  depth  resolved  cathodolumlnescence  Is  discussed  and  the 
proper  action  taken  to  minimize  its  impact  is  given.  The  L(V)  curves  are 
presented  in  the  third  section.  In  the  final  section  the  analytical  and 
experimental  L(V)  curves  are  compared.  The  physical  parameters  chosen 
for  the  analytical  curves  are  justified  based  on  other  workers'  results. 

Conclusions  based  on  this  work  and  recommendations  for  future  work 
are  presented  in  the  last  chapter.  The  main  conclusion  is  that  a  basis 
for  quantitative  catliodoluminescence  has  been  established.  The  recommen¬ 
dations  generally  concern  Improvements  in  the  experimental  facilities  and 
the  theoretical  model. 


LI.  The  Model 


In  order  to  analytically  determine  the  luminescent  Intensity  as  a 
function  of  the  electron  beam  energy,  the  following  processes  need  to  be 
considered.  The  electron  beam  loses  energy  to  the  GaAs  by  creating  elec¬ 
tron-hole  pairs;  thus  the  method  of  calculating  the  rate  of  energy  loss 
with  depth  is  important  here  as  this  determines  the  generation  function 
for  electron-hole  pairs.  The  electron-hole  pairs  just  created  then 
diffuse  through  the  crystal  before  recombining.  As  far  as  the  lumines¬ 
cence  curve  is  concerned,  the  position  where  the  electrons  and  holes 
recomoine  is  important,  not  the  position  where  they  were  generated.  The 
actual  luminescence  that  emerges  from  the  surface  can  then  be  calculated. 
This  requires  that  both  absorption  of  the  recombination  radiation  and  the 
density  of  recombination  centers  be  considered. 

The  luminescence  versus  beam  voltage,  L(V),  curve  can  be  calculated 

by  either  a  completely  analytical  approach,  or  by  a  numerical  approach. 

( 18} 

The  analytical  approach  has  been  followed  by  Gergely'  ',  Wittry  and 
Kyser^’  Rao  Sahib  and  Wittry^^,  and  Fano^^.  Gergely  assumed 

a  depth-dose  curve  and  solved  the  diffusion  equation  for  the  excess 
carrier  density.  Given  that  the  recombination  centers  were  uniform  in 
depth,  he  was  able  to  Integrate  the  excess  carrier  expression  to  find  the 
L(V)  curve. 

Wlttry/Kyser  and  Rao  Sahib/Wittry  determined  L(V)  curves  by  solving 

the  diffusion  equation  using  the  Green's  function  determined  by  van  Roos- 

(22) 

broeck  .  In  reference  (19)  the  luminescence  intensity  is  assumed  pro¬ 
portional  to  the  excess  carrier  density  raised  to  some  power  between  1 
and  2,  while  in  reference  (9)  the  power  is  assumed  to  be  1. 
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(21) 

Fano  takas  a  different  approach.  He  solves  for  the  integral  of 

the  product  of  the  energy  lost  by  electrons  within  the  crystal  and  the 
average  depth  of  the  energy  loss.  He  assumes  a  model  in  which  surface 
effects  are  more  important  than  bulk  effects  and  comes  up  with  a  Boltz- 
mann-llke  equation  that  he  can  solve  analytically. 

The  numerical  solutions  are  based  on  either  a  numerical  solution  of 
the  Boltzmann  equation^'*’  or  a  Monte  Carlo  approach^"*’ 

None  of  the  authors  who  used  numerical  techniques  calculated  L(V)  curves. 
As  they  were  all  motivated  by  electron  microprobe  work,  they  usually 
calculated  the  spatial  distribution  of  x-ray  production  in  solid  tar¬ 
gets.  Their  principal  result  of  interest  is  the  calculation  of  the  rate 


at  which  electrons  lose  energy  as  they  penetrate  a  solid  target  (called 


energy  loss  or  depth-dose  curves). 


The  Boltzmann  equation  approach  takes  the  transport  equation. 


(24,  39) 


(sinOf (x,0 ,s) ) 


cost) i — [s  l  n*3  f  ( x ,  6  ,  s )  ] 

a  X 


i _  3 

X(s)  30* 


in&gglf (x,0,s)] 


and  so  Ives  tor  f(x,t),a),  the  electron  distribution  function.  In  this 
equation,  x  Is  the  depth  Into  (he  crystal,  s  is  the  path  length  of  the 
electron,  0,  the  angle  the  electron  makes  with  the  normal  to  the  surface, 
and  X  is  the  transport  mean  free  path.  A  complete  solution  for  f(x,0,s) 
gives  the  complete  history  of  electron  transport  in  the  specimen.  Since 
energy  depends  only  on  the  path  length  traveled  and  not  on  the  depth, 
the  path  length,  s,  gives  I  he  energy  ot  the  beam  at  depth  x  headed  In 
the  U  direction.  From  here  the  depth-dose  curve  can  be  calculated,  then 
the  diffusion  equation  solved  for  the  excess  carrier  density  and  the 
L(V)  curve  calculated. 
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The  Monte  Carlo  approach  uses  a  Monte  Carlo  technique  to  determine 
the  depth-dose  function.  Using  the  depth-dose  curve  as  a  forcing  func¬ 
tion,  the  diffusion  equation  is  solved  for  the  excess  carrier  density. 
The  L(V)  curve  can  be  found  by  integrating  the  product  of  the  excess 
carrier  density  and  the  recombination  center  density  over  depth.  It  is 
possible  to  solve  the  problem  through  to  the  L(V)  curve  using  Monte 
Carlo  techniques. 

An  analytical  approach  was  rejected  for  this  study  since  the  para¬ 
meters  on  electron  beam  penetration  into  GaAs  are  not  readily  available. 
This  means  that  either  a  Monte  Carlo  simulation  or  a  transport  equation 
solution  would  be  required  to  provide  the  necessary  information  to  eval¬ 
uate  the  theory.  The  Monte  Carlo  equation  approach  was  selected  because 
there  exists  a  bigger  body  of  literature  and  the  results  are  better 
documented  than  for  the  Boltzmann  equation  approach.  The  Monte  Carlo 
analysis  also  has  the  advantage  of  being  more  conceptually  transparent 
than  the  transport  equation  approach.  In  addition,  this  approach  makes 
data  interpretation  and  program  modification  easier.  Of  the  two  Monte 
Carlo  approaches,  i.e.  total  or  partial  simulation,  the  total  Monte 
Carlo  simulation  was  ruled  out  because  the  scattering  cross  sections 
necessary  are  not  readily  available  and  because  the  partial  simulation 
requires  much  less  computer  time  to  execute. 

This  chapter  is  divided  into  three  parts.  In  the  first  part,  the 
Monte  Carlo  analysis  of  electrons  penetrating  into  GaAs  is  developed. 
This  includes  a  section  describing  the  past  work  on  electron  penetration 
into  solids,  a  section  describing  how  the  Monte  Carlo  calculation  was 
performed,  and  a  section  of  results  where  the  Monte  Carlo  analysis  is 


compared  to  other  workers'  results.  In  the  next  part  of  the  chapter, 
the  diffusion  equation  Is  solved  for  the  excess  carrier  density  and  the 
L(V)  curve  is  calculated.  The  appropriate  diffusion  equation  is  first 
derived  and  then  the  computer  programs  that  are  used  to  evaluate  the 
excess  carrier  density  and  the  L(V)  curves  are  described.  Finally,  the 
results  of  the  L(V)  calculation  are  presented  in  the  form  of  graphs  that 
show  how  the  L(V)  curve  varies  as  the  diffusion  length,  surface  recom¬ 
bination  velocity,  absorption  coefficient,  and  impurity  profile  are 
changed.  In  the  last  part  of  this  chapter  the  calculation  procedure  is 
summarized  for  future  easy  reference. 

Electron  Penetration  into  GaAs:  Background 

The  theoretical  analysis  of  electron  beam/target  interaction  falls 
into  three  categories.  These  are:  single,  plural,  and  multiple  scat¬ 
tering  of  the  electron  beam,  depending  on  whether  there  are  one,  a  few 
(less  than  25),  or  many  electron-atom  interactions.  This  usage  differs 
from  that  of  high  energy  scattering  where  single,  plural,  and  multiple 
scattering  are  synonymous  with  Large,  medium,  and  small  angle  scattering 
respectively. 

Single  scattering  results  from  elastic  collisions  between  electrons 
and  atomic  nuclei.  The  theoretical  treatment  starts  with  the  classical 
Rutherford  formula.  Various  theories  exist  to  take  account  of  the  screen¬ 
ing  effect  of  the  electron  cloud  for  the  nucleus.  The  most  frequently 
applied  corrections  to  Rutherford  scattering  are  the  Wentzel  and  the 
Thoroas-Fermi  statistical  model.  The  Wentzel  model  assumes  an  exponential 
fall  of  the  screening  field.  The  Thomas-Fermi  approach  determines  the 
field  for  the  atoms  by  considering  the  atomic  electrons  as  a  degenerate 
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gas.  The  potential  so  obtained  is  used  within  the  first  Born  approxima¬ 
tion  to  calculate  the  differential  scattering  cross  section.  Everhart 
used  the  single  scattering  assumption  to  develop  a  simplified  theory  for 
the  reflection  of  electrons  from  solids.  The  theory  of  plural  scattering 
has  been  discussed  by  Bothe^1,  ^  ,  Wentzel^^,  Moliere^\  Lenz^35\ 
and  Smith  and  Burge^^  among  others.  CossJett  and  Thomas  have 
reviewed  these  treatments  in  an  attempt  to  explain  their  experimental 
results.  The  experimental  work  reported  by  them  in  reference  37  showed 
good  agreement  with  Bothe's  theory  in  the  form  due  to  Lenz.  Bothe  had 
been  able  to  give  a  formal  solution  for  plural  scattering  on  the  basis 
of  a  statistical  approach  called  error  theory.  In  Bothe's  analysis  he 
assumed  successive  collisions  to  be  statistically  independent  and  ne¬ 
glected  large  angle  single  scattering  and  energy  loss.  The  problem  with 
Bothe's  theory  is  that  the  two  main  equations  are  very  difficult  to 
integrate.  Lenz  was  able  to  integrate  the  first  of  Bothe's  equations 
and  cast  the  second  one  in  a  form  for  numerical  integration. 

(32 

The  multiple  scattering  theories  were  first  developed  by  Bothe  ’ 

38)  (39) 

,  Bethe,  Rose  and  Smith  and  on  slightly  different  lines  by  Goudsmit 

and  Saunderson^^  ,  Lewis ,  Spencer^^’  and  Meister^^.  Cosslett 
(45) 

and  Thomas  have  also  discussed  the  multiple  scattering  theories, 
where  they  compared  the  various  approaches  against  measurements  on 
aluminum,  copper,  silver  and  gold.  Bothe's  approach  is  based  on  error 
theory  as  is  his  plural  scattering  theory.  The  others  are  based  on  the 
diffusion  equation.  Bethe,  Rose,  and  Smith  neglect  energy  loss  and  use 
the  Fokker-Planck  approximation  to  the  diffusion  equation  which  neglects 
large  angle  single  scattering.  Goudsmit/Saunderson,  Lewis,  Spencer,  and 
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Meister  included  the  energy  loss,  but  only  considered  infinite  targets 

where  the  source  was  assumed  embedded  in  the  target.  These  theories 

were  moderately  successful  for  high  beam  energies  where  the  common 

assumption  of  small  angle  scattering  holds.  The  attempts  to  modify  the 

theory  for  energies  in  the  1  to  50  keV  range  were  by  Moliere^^  and 
(47) 

Lenz  .  All  of  these  theories  assume  that  the  beam  is  normally  inci¬ 
dent  to  the  target  surface. 

With  the  advent  of  larger  computers,  numerical  solutions  to  the 

penetration  problem  have  been  more  prevalent  lately.  Brown  and  Ogil- 

vie^^  solve  the  Boltzmann  transport  equation  numerically  using  the 

multiple  scattering  approach  of  Bethe,  Rose  and  Smith.  Later  Brown, 

(49) 

Wittry  and  Kyser  refined  the  calculations  by  considering  different 

approximations  over  different  path  legnths.  For  instance,  they  use  a 
single  scattering  approximation  to  the  Boltzmann  transport  equation  for 
short  path  lengths.  The  result  of  their  calculation  is  a  distribution 
function  f(x,0,S)  that  gives  the  probability  that  an  electron  which  has 
traveled  a  distance  between  S  and  S  +  dS  will  have  a  position  between  x 
and  x  +  dx  below  the  specimen  surface  while  at  the  same  time  its  direc¬ 
tion  of  travel  makes  an  angle  with  the  internal  normal  to  the  specimen 
surface  of  between  0  and  0  +  dO.  The  distribution  function  is  then  used 
to  calculate  physically  observed  results.  An  alternate  approach  is  a 
Monte  Carlo  simulation.  There  are  two  kinds  of  simulations  that  may  be 
employed.  Direct  simulation  considers  each  collision  an  electron  would 
have  as  it  penetrates  the  target.  The  alternative  is  a  "condensed" 
random  walk  where  many  electron- target  collisions  are  condensed  into  one 
equivalent  collision.  Calculations  of  the  latter  have  been  performed 
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by  Green^^,  Bishop^^  and  Shimizu  et  al^^.  Direct  simulations  have 
been  made  by  Shimizu  et  al^^  and  by  Green  and  Leckey^^.  A  detailed 
explanation  of  the  interaction  of  1  to  50  keV  electrons  with  solids  has 
yet  to  be  achieved  despite  all  of  the  theoretical  activity  cited  above. 
One  of  the  limitations  in  the  application  of  the  theoretical  models  to 
scattering  experiments  from  solids  lies  in  their  reliance  on  adjustable 
parameters  and  in  the  fact  that  they  provide  little  insight  into  the 
dynamics  of  the  individual  collision  process  involved.  A  purely  analy¬ 
tical  approach  to  the  problem  of  electron  scattering  based  on  a  consid¬ 
eration  of  individual  scattering  events  would  present  formidable  pro¬ 
blems.  Realistic  boundary  conditions  are  also  very  difficult  (if  not 
impossible)  to  handle  theoretically.  On  the  other  hand,  these  problems 
are  easily  handled  with  a  Monte  Carlo  analysis.  The  best  agreement 
between  experiment  and  theory  for  electrons  penetrating  into  solids 
comes  from  a  Monte  Carlo  analysis  (see  references  53  and  54).  In  addi¬ 
tion,  a  Monte  Carlo  analysis  has  the  Intuitive  appeal  of  being  just  as 
the  theoretician  pictures  the  experiment  occurring.  Thus,  the  results 
tend  to  be  more  transparent  than  In  more  obscure  theoretical  approaches. 

Although  a  great  deal  of  effort  has  been  put  into  studying  elec¬ 
tron-beam  penetration  into  solids  comparitively  little  work  has  been 
done  with  crystalline  solids.  References  that  apply  to  GaAs  are  Stim- 
ler^^,  Anderson^^,  Kiein^^  ,  Schiller  and  Boulou^"^  ,  Norris  et 
ai^13\  and  Wittry  and  Kyser^60\  Slimier  gives  families  of  curves  of 
the  depth  of  electron  penetration  versus  density  of  the  solid  with 
electron  energy  as  a  parameter.  These  curves  are  deduced  from  the  data 
of  Ehrenberg  and  King^"^.  He  only  considers  electron  beams  normal  to 
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the  surface  and  does  not  plot  "depLh-dose"  curves  (energy  dissipation 
versus  depth  Into  the  target).  Anderson  applies  Spencer's  '  theore¬ 
tical  work  for  electrons  released  in  an  infinite  medium  to  GaAs.  While 
the  boundary  conditions  are  different  (the  target-vacuum  interface  in 
Anderson's  case  is  at  depth  equals  zero  while  Spencer's  theory  is  for  an 
infinite  medium),  Anderson  calculates  depth  dose  curves  for  25,  50  and 
100  keV  beams  anyway.  He  does  not  consider  the  effect  of  different 
angles  of  incidence  on  the  target  nor  does  he  indicate  what  to  do  for 
voltages  below  25  keV  (Spencer's  lowest  energy). 

Klein  has  attempted  to  develop  a  phenomenological  model  capable  of 
ciescribing  all  pertinent  aspects  of  electron-beam  penetration  that  will 
agree  with  experimental  evidence  wherever  comparisons  can  be  made. 

Klein  assumes  the  electrons  penetrate  into  the  crystal  without  loss  to  a 
depth,  R^,  the  depth  of  complete  diffusion,  where  they  diffuse  randomly 
in  all  directions,  transferring  energy  to  the  medium  at  an  exponentially 
decreasing  rate  along  any  radius  vector  of  the  sphere  of  excitation  (see 
fig.  1).  Klein's  phenomenological  curves  reproduce  Anderson's  calcula¬ 
tions!  results  for  GaAs.  Schiller  and  Boulou  use  Klein's  model  to 
calculate  depth-dose  curves  at  beam  energies  lower  than  25  keV.  Again 
Klein's  data  Ignores  the  effect  of  the  angle  of  electron  beam  incidence 
on  the  target.  Also  ignored  is  the  crystalline  structure  of  the  target 
since  the  results  are  the  same  for  amorphous  as  well  as  crystalline 
materials . 

(13) 

Norris  et  al  present  depth  dose  curves  for  electrons  penetrat¬ 
ing  Into  GaAs  at  45°  incidence  wilh  energies  of  5,  10,  and  20  keV. 

These  curves  were  obtained  by  scaling  an  experimentally  determined  high 
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energy  (2  MeV)  curve  to  the  low  energy  (a  few  keV)  normal  incidence 
electron  range.  Their  justification  for  this  procedure  is  that  the 
depth-dose  curve  is  approximately  a  universal  curve  that  can  be  scaled 
given  a  known  electron  range  for  a  given  incident  electron  energy.  This 
has  only  been  established  at  keV  energies  and  normal  incidence  but  they 
assume  that  the  high  energy  curves  can  be  scaled  to  the  keV  range  for 
various  angles  of  incidence.  Wittry  and  Kyser  assume  that  the  depth- 
dose  curve  is  Gaussian.  They  consider  only  normal  incidence  at  29  keV 
and  determine  the  mean  and  standard  deviation  from  the  transport  cal¬ 
culation  made  for  them  by  Brown  following  the  technique  in  reference 
(48) .  They  suggest  that  the  parameters  for  other  energies  be  obtained 
by  scaling  the  29  keV  data,  by  taking  experimental  data,  or  by  calcula¬ 
ting  using  the  transport  equation  approach  or  a  Monte  Carlo  technique. 
The  Monte  Carlo  technique  is  the  approach  chosen  here. 

Monte  Carlo  Analysis 

This  section  will  begin  with  a  summary  of  the  Monte  Carlo  calcula¬ 
tion,  with  certain  parts  of  the  calculation  developed  in  detail.  Then 
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a  brief  description  of  the  computer  programs  used  to  perform  the  analy¬ 
sis  will  be  given. 

The  electron  energy  loss  curves  (called  depth-dose  curves)  for 
electron  beams  penetrating  into  GaAs  were  calculated  using  the  standard 
Monte  Carlo  procedures  as  reviewed  by  Berger  This  procedure  con¬ 

sists  of  condensing  several  elastic  atomic  collisions  into  one  equiva¬ 
lent  collision.  The  scattering  angle  is  randomly  chosen  from  a  table  of 
angles  that  is  constructed  from  the  multiple  scattering  theory  of  Goudsmit 
and  Saunderson^*^  .  The  length  or  amount  of  material  the  beam  must  go 
through  to  suffer  the  required  number  of  elastic  collisions  is  called 
the  step  size.  The  particular  formulation  used  here  takes  into  account 
the  step  size  and  the  energy  of  the  beam  aL  the  beginning  of  each  step. 
Bethe's  law  as  given  in  reference  (62)  is  used  to  account  for  the  energy 

loss  in  each  step.  The  step  size  and  location  of  the  point  of  the 

(53) 

scattering  follow  the  method  of  Shimizu  et  al  .  In  this  method,  the 
step  size  is  scaled  according  to  the  energy  of  the  beam  at  the  beginning 
of  each  step,  the  initial  beam  energy,  and  the  initial  step  size.  Given 
the  angle  through  which  the  electron  is  going  to  scatter,  the  location 
of  the  scattering  point,  the  energy  of  the  electron  at  the  beginning  of 
the  step,  and  the  energy  lost  in  the  step,  the  trajectory  and  energy 
loss  over  that  trajectory  can  tie  calculated.  The  energy  loss  curve  is 
then  calculated  by  dividing  the  thickness  of  the  target  into  equal 
increments,  called  bins,  and  then  calculating  how  much  energy  is  lost  in 
each  bin.  At  this  point,  a  mure  detailed  discussion  will  be  given  on 
the  angular  distribution  function,  the  step  size,  the  method  of  random 
sampling  and  Bethe's  energy  loss  law. 
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Goudsmit  and  Saunderson^*^  studied  the  angular  distribution  of 
multiple-scattered  particles.  They  derived  the  exact  angular  distribu¬ 
tion  function  as  a  Legendre  series, 

A(co)sino)do)  *=  £(k  +  -2)expl-/oG^(s  )ds  }P^(cosa3)sin(Ddco  (1) 

k=o 

where  G  (s)  =  2ttN/\i  (0  ,  s)  { 1  -  P,  (cos8) } sin0d0 . 

K.  O  K. 

N  is  the  number  of  atoms  per  unit  volume,  s  is  the  path  length  traveled 
by  the  particle  in  that  particular  step,  and  o(0,s)  is  the  single  scat¬ 
tering  cross  section.  Since  the  electrons  do  not  lose  much  energy  over 
any  one  step,  the  collision  cross  section  is  assumed  constant  over  that 
particular  step.  Thus  in  equation  (1), 

fSG,  (s')ds'  =  /ASG.  (s')ds'  «  C.  (s) As 
ok  ok  k 

The  Goudsmit-Saunderson  distribution  applies  to  all  angular  deflec¬ 
tions  regardless  of  their  magnitude.  It  can  also  be  evaluated  for  any 
desired  single  scattering  cross  section.  The  cross  section  that  is  used 
here  is  the  screened  Rutherford  cross  section 

a(0,s)  -  (Z2eA)/(p2v2(l  -  cosO  +  2n)2)  (2) 

where  Z  is  the  equivalent  atomic  number  (32  for  GaAa) ,  e  is  the  electron 
charge,  p  is  the  electron  momentum,  v  is  the  electron  velocity,  and  r|  is 
a  parameter  that  takes  into  account  the  screening  of  the  nuclear  charge 
by  the  orbital  electrons,  n  is  considered  an  adjustable  parameter  whose 
value  can  be  approximately  determined  from  a  formula  by  Nigam  et  al^^, 

n  »  !itl.l2(h/p)(Z1/3/0.88ba,)}2  (3) 

where  a.  is  the  Bohr  hydrogen  radius  and  k  is  Planck's  constant. 
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The  path  length,  s,  that  the  electron  will  travel  Is  broken  down 


into  a  number  of  steps,  AS^.  The  step  size  over  which  the  angular  dis¬ 
tribution  Is  to  be  calculated  must  be  long  enough  to  satisfy  the  treat¬ 
ment  of  multiple  scattering  and  to  obtain  stable  convergence  of  equation 
(1).  On  the  other  hand,  AS^  should  be  as  short  as  possible  to  improve 
the  accuracy  of  the  Monte  Carlo  calculation.  Shinoda  et  al^"^  have 
shown  that  by  choosing  the  i*"*1  step  size  to  be 


AS.  =  -=r  AS 
i  E  o 
o 


(A) 


,  th 


where  E^,  the  beam  energy  at  the  beginning  of  the  i  step,  is  given 
by 


P  —  F  4- 

i  "  i-1  ds 


AS. 
i 

E  =  E 


(5) 


i-1 


the  accuracy  of  the  calculation  is  about  the  same  in  each  step.  The 
quantity  ASq  is  the  first  step  length  and  must  be  given.  This  initial 
step  size  is  hard  to  determine.  Its  value  can  be  estimated  from  Coss- 
lett  and  Thomas  or  taken  as  approximately  the  value  for  copper  found 

/  C  O  \ 

by  Shinoda  et  al'  .  As  the  initial  step  size  is  a  parameter  that  is 
to  be  adjusted  to  improve  agreement  with  experimental  results,  Shinoda's 
value  for  copper,  .18  microns  at  30  keV  beam  energy,  will  be  used. 

The  procedure  used  to  sample  random  variates  from  a  distribution 
function  with  the  use  of  pseudorandom  numbers  relies  on  the  calculation 
of  the  cumulative  probability  distribution  Junction  (see  any  book  on 
Monte  Carlo,  e.g.  references  50  or  66) 


where  f(x)  is  the  distribution  function  to  be  sampled.  In  order  to 
sample  f(x): 


(a)  Compute  F(x)  for  a  dense  set  of  x  values  over  the  range  of 


f  (x)  . 


(b)  By  interpolation,  find  a  set  of  x  values  such  that 


F(X  )  -  (m  -  lj)/M  m  =  1,2.  .  .M 

m 

(c)  Store  the  F  ^(X  )  in  computer  memory. 

m 


(6-A) 


(d)  Choose  a  random  number  p. 

(e)  The  desired  f(x)  is  the  one  located  at  memory  location  i  = 
integral  part  of  pM. 

The  Coudsmit  and  Saunderson  distribution  is  used  to  determine  the  polar 
angle.  The  azimuthal  angle  is  distributed  from  0  to  2 u  provided  the 
medium  is  isotropic  and  polarization  is  ignored. 

The  electrons  are  assumed  to  lose  energy  continuously.  This  is 
called  the  continuous  slowing  down  approximation.  In  this  approxima¬ 
tion,  the  effects  of  straggling  are  ignored.  Energy  loss  is  usually 


given  in  terms  of  Bethe's  law 


(67) 


The  particular  form  chosen  for 
(62) 


Bethe's  law  is  given  by  Berger  and  Seltzer''  '  following  the  formulation 
of  Rohrlich  and  Carlson^*^, 


dE 

ds 


2ttn  pr2mc2  „ 
3  z 


3 


-  log  T-  ^T-  +22\  +  F-(t)  -  5 
A  2(l/mc  ; 


(7) 


F-(t)  =  1  -  B2  +  [t2/8  -  (2t  +  l)log2]/(x  +  l)2 


where 


me 


rest  energy  of  electron  ■  0.511  MeV 
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X  =  kinetic  energy  in  units  of  me 
3  =  [ v ( i  +  2)]  V (t  +  1)  =  velocity /c 

Z  =  atomic  number 
A  =  atomic  weight 
p  =  density 

I  =  mean  excitation  energy 
5  =  density  effect  correction 


=  Avagardro's  number 


r2  =  ,  2,  2,2 


r‘  =  (e'/mcT  =  7.904030  x  10  26  cm2, 
o 


The  density  effect  correction,  6,  takes  into  account  the  reduction  of 
the  collision  loss  due  to  polarization  of  the  medium.  Based  on  the 
values  given  for  copper,  the  density  effect  correction  is  negligible  at 
these  beam  energies.  The  mean  excitation  energy,  I,  is  chosen  to  con¬ 


form  to  the  recommendations  of  Berger  and  Seltzer 


(62) 


They  give  an 


approximate  formula  for  I  which  they  call  1^  since  it  is  slightly 

larger  than  the  I  defined  in  terms  of  oscillator  strengths.  It  is  the 

I  ,,  value  that  is  used  for  this  Monte  Carlo  calculation, 
adj 


I  ,  =  Z  (9.76  +  58.8  Z-1'19) 

adj 


(8) 


For  mixtures  and  compounds,  the  mean  energy  loss  is  assumed  to  be  the 
sum  of  the  losses  in  the  constituent  elements.  Thus 


and 


-1  Z 

108  Vi  "  f  I5rj-pil08ladj,j 


z/A  ■  ?  5  va 


(9) 


For  GaAs,  Z  =  32,  p  =  5.32  gr.  and  I  ,,  =  342.4  ev. 

adj 
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Civen  the  above  values,  the  cumulative  distribution  function  for 
multiple  scattering  of  electrons  in  GaAs  is  calculated.  The  program 
that  does  this  is  called  LMSD.  It  is  used  to  evaluate  a  cumulative  dis¬ 
tribution  function  for  each  as  determined  from  equation  (5)  down  to 
an  E^  approximately  equal  to  1  keV.  At  this  point,  the  velocity  of  the 
incoming  electron  becomes  about  the  same  as  the  velocity  of  the  atomic 
electrons.  The  electron  stopping  formula  [equation  (7)]  is  no  longer 
valid.  Since  the  exact  form  is  not  known,  an  arbitrary  interpolation 
procedure  is  used.  Berger  and  Seltzer  follow  Nelms^^  and  assume 
(dE/ds)  =  0  at  E^  =  0  and  linearly  interpolate  from  E^  =  0  to  =  1 
keV.  Here  the  same  procedure  is  used.  Next,  a  table  of  angles  is  con¬ 
structed  according  to  formula  (6-A).  These  are  stored  on  data  cards  for 
use  in  the  next  phase  of  the  calculation. 

The  program  used  to  calculate  the  energy  loss  versus  depth  curves 
is  called  EPIGA  for  electron  penetration  into  GaAs.  The  main  parts  of 
the  program  are  shown  in  figure  (2) .  The  program  uses  the  cumulative 
distribution  functions  from  LMSD  and  the  GaAs  parameters  as  inputs  and 
calculates  the  energy  loss  versus  depth  curve  called  the  DEDS  curve.  Its 
output  la  a  table  of  the  energy  an  average  electron  is  expected  to  lose 
at  that  depth  versus  depth.  Included  in  the  output  is  the  average 
energy  loss  versus  depth  for  the  backscattered  as  well  as  absorbed 
beams.  The  number  of  electrons  backscattered  and  absorbed,  the  total 
number  of  electrons,  the  angle  of  Incidence  for  the  beam,  the  incident 
energy  and  total  beam  energy  are  then  compared  to  make  sure  all  of  the 
energy  is  accounted  for. 
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The  intermediate  steps  in  the  e a  1 eu 1  a L i un  inelude  cumulating  t he 
energy  loss  versus  path  length,  the  electron  trajectory  and  the  amount 
ot  energy  lost  in  each  bin.  1  he  energy  loss  versus  path  length  is  the 
same  for  all  electrons.  This  is  just  another  way  of  saying  that  strag¬ 
gling  has  been  ignored.  This  does  not  mean  that  all  electrons  have  the 
same  energy  loss  at  a  particular  depth.  What  it  does  mean  is  that  all 
electrons  that  have  traveled  through  the  same  amount  of  GaAs  have  lost 
the  same  amount  of  energy  regardless  of  their  respective  depths.  Next, 
a  matrix  Is  formed  that  lines  up  the  depth  in  the  material  with  the 
depth  at  which  the  particle  scattered  and  the  energy  loss  over  that 
step.  The  axis  normal  to  the  target  surface  is  divided  up  into  equal 
increments  called  bins.  Subroutine  Score  is  then  used  to  take  this 


Fig.  2  Calculation  of  the  Energy  Loss  Curve 
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matrix  and  decide  how  much  energy  is  lost  in  each  bin.  A  number  of 
electron  trajectories  are  computed  and  the  average  energy  loss  versus 
depth  curve  is  calculated. 

Results  from  the  Program 

In  order  to  verify  the  EPIGA  computer  program,  the  results  for 
electrons  at  30  keV  normally  incident  on  Cu  were  checked  against  those 
of  Shimizu^^.  The  results  of  this  check^^  verified  that  EPIGA  is 
functioning  properly. 

As  mentioned  earlier  n  and  ASq  are  considered  adjustable  parameters 

that  help  match  the  results  from  the  model  to  experimental  results. 

Unfortunately,  there  is  not  a  great  deal  of  information  on  electron  beam 

interaction  with  GaAs  as  contrasted  to  the  work  on  Al,  Cu,  Ag,  and  Au. 

One  data  point  that  is  available  is  the  one  reported  by  WitLry  and 

Kyser^^.  They  used  a  backscattered  fraction  of  0.33  for  a  30  keV  beam 

of  electrons  normally  incident  on  GaAs.  A  number  of  runs  with  different 

combinations  of  n  and  AS  showed  that  AS  *  .185  micron  and  n  ■  .0054972 

o  o 

reproduce  Wittry  and  Kyser's  result.  The  Monte  Carlo  simulation  pre¬ 
dicted  a  backscattered  fraction  of  33.78%.  An  additional  checkpoint  is 
given  by  Klein^^  .  He  calculates  the  fraction  of  beam  energy  that 
should  be  backscattered  for  various  elements  and  compounds.  For  a  30 
keV  beam  of  electrons  normally  incident  on  GaAs,  he  gets  an  energy  loss 
fraction  of  23.7%.  The  Monte  Carlo  result  is  22.9%  which  agrees  very 
well. 

Figure  (3)  shows  depth  dose  curves  for  a  beam  incident  at  45°  with 
energies  of  5,  10,  15,  20  keV.  Figure  (4)  is  a  least  squares  approxima¬ 
tion  by  cubic  splines  with  variable  knots.  The  curve  fit  Is  performed 

24 


5 


by  subroutine  ICSVKV,  an  1MSL  subroutine.  The  curve  is  truncated  when¬ 
ever  the  tail  of  the  least  squares  fit  goes  negative  or  turns  up.  The 
curve  is  also  terminated  when  the  raw  data  runs  out.  These  curves  can 

(13) 

be  compared  to  energy  loss  curves  for  Norris  et  al  shown  in  figure 


0  0  0.2  0.4  0 .6  OS  10  li  14  It 
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(13) 

Fig.  5  Depth-Dose  Curves  -  Norris  et  al 

The  Monte  Carlo  analysis  predicts  a  deeper  penetration  than  the  curves 
Norris  derives.  Martinelli  and  Wang^^  measured  the  depth  that  elec¬ 
tron  beams  of  normal  incidence  penetrate  into  GaAs.  In  their  experi¬ 
ment,  GaAs  thin  films  were  grown  on  MgAl^O^  spinel  substrates.  When  the 
electrons  had  sufficient  energy  they  caused  the  MgAl^O^  to  glow  when 
they  struck  it.  liy  noting  at  what  voltage  the  electron  beam  just  pene¬ 
trated  through  the  GaAs  Into  the  MgAl,;0^,  and  by  knowing  the  thickness 
of  the  thin  film,  as  well  as  correcting  for  the  residual  energy  left  in 
the  beam  in  order  to  make  the  MgAl^O^  glow,  they  were  able  to  establish 


a  depth  of  maximum  penetration.  Table  I  compares  their  results  witli 


those  of  this  study.  In  all  cases  their  range  is  greater  than  the  one 
predicted  here.  There  are  several  reasons  for  this.  One  reason  is  that 
the  electrons  had  to  have  some  energy  left  when  they  struck  the  MgAl^O^ 
substrate  in  order  to  be  detected.  While  Martinelli  and  Wang  corrected 
their  data  for  this,  the  lower  energy  scattering  is  more  isotropic  than 
the  higher  energy  scattering,  the  correction  Martinelli  and  Wang  apply 
overcorrects  their  range  values,  making  them  longer  than  they  should  be. 
Another  reason  is  that  the  Monte  Carlo  data  at  the  deepest  penetration 
depths  is  subject  to  large  statistical  errors;  there  are  not  many  elec¬ 
trons  that  make  it  to  this  depth.  Finally,  the  lower  energy  loss  rate 
was  arbitrarily  set.  While  this  does  not  change  the  overall  shape  of 
the  depth-dose  curve  much,  it  could  affect  the  tail  of  the  curve.  In 
spite  of  these  problems,  the  agreement  is  considered  excellent. 

Table  I.  Comparison  of  Martinelli  and  Wang's  Experimental 
Penetration  Data  to  the  Monte  Carlo  Analysis 


Energy  (keV) 

M  &  W  Maximum 
Depth  (microns) 

Monte  Carlo  Maximum 
Depth  (microns) 

%  Difference 

5 

.283 

.250 

11.6% 

10 

.779 

.725 

6.9% 

15 

1.-408 

1.325 

5.9% 

20 

2.142 

2.125 

.8% 

The  minor  disagreement  with  Martinelli  and  Wang's  data  does  not  affect 
the  rest  of  the  calculation  very  much  since  at  this  depth  the  energy 
loss  curve  lias  dropped  to  a  point  where  very  few  electron/hole  pairs  are 
being  generated  compared  to  shal lower  depths.  Figure  (6)  shows  the  depth- 
dose  curves  for  5,  iO,  20  and  30  keV  beams  at  normal  incidence. 


Finally,  figures  (7),  (8)  and  (9)  are  of  the  energy  loss  curves  at 


45°  incidence  for  the  remaining  beam  voltages  used  in  the  diffusion 
calculation. 

Calculation  of  the  Excess  Carrier  Density  and  Recombination  Radiation 

The  electron  beam  gives  up  its  energy  to  the  GaAs  crystal  by  creating 
electron/hole  pairs  and  phonons.  The  electron/hole  pairs  then  diffuse 
away  from  the  creation  point,  recombining  at  some  other  place  in  the 
crystal.  Some  of  the  recombination  radiation  is  given  up  as  luminescence. 
The  goal  of  this  study  is  to  be  able  to  calculate  how  the  luminescent 
intensity  will  vary  with  beam  energy.  if  the  excess  carrier  density, 
density  of  available  recombination  states,  and  some  exponential  absorption 
factor  for  the  emitted  radiation  are  known,  the  L(V)  curve  can  be  calcu¬ 
lated.  How  that  is  done  will  be  discussed  in  this  section.  The  first 
part  is  a  derivation  of  the  appropriate  diffusion  and  L(V)  equation. 

Next  comes  an  outline  of  the  computer  programs  that  were  written  to 
solve  these  equations.  This  is  followed  by  a  results  section  that  shows 
how  the  L(V)  curves  change  with  various  parameter  changes.  Last  of  all 
is  a  summary  of  the  L(V)  calculation,  reviewing  briefly  how  all  the 
parts  fit  together. 

Justification  of  Approach 

The  approach  taken  is  to  so Lve  the  one-dimensional  diffusion  equation 
for  the  excess  carrier  density  and  then  solve  for  the  L(V)  curve.  In 
this  section,  the  approximations  that  are  made  in  doing  this,  how  likely 
they  are  to  be  met,  and  in  some  cases,  what  alternatives  are  available, 
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Any  elementary  book  on  solid  state  theory  or  semiconductor  physics 

(72) 

will  contain  a  discussion  of  the  diffusion  equation  (see  McKelvey 

(79) 

for  instance).  The  diffusion  equation  can  be  written  as  ' 


-V  •  J  +  g  -p/x  =6p/6t  (10) 

P  P  P 


-V  •  J  +g  -n/x  =  6n/6t 

(ID 

n  n  n 

with 

J  =  -D  Vp  +  py  E 

(12) 

P  P  P 

J  =  -D  Vn  -  np  E 
n  n  n 

(13) 

the  subscripts  p  and  n  refer  to  holes  and  electrons  respectively.  The 
J’s  are  particle  flux  densities,  the  g's  are  generation  rates,  and  the 
x's  are  recombination  rates.  Equations  (10)  and  (11)  are  frequently 
called  the  continuity  equations.  When  (12)  and  (13)  are  substituted 
into  (10)  and  (11),  the  results  are 


D  V2p  +  Vp 
P 

D  V2n  +  Vn 
n 


VD  -  V  *  (p  pE)  +  g  -  p/x  =*  6p/6t  (14) 

P  P  P  P 

VD  +  V  •  (p  nE)  +  g  -  n/x  ■  6n/6t  (15) 

n  n  n  n 


In  equations  (14)  and  (15)  there  is  a  term  which  depends  on  the  depth 

dependence  of  the  diffusion  constant.  In  general,  a  variation  in  the 

diffusion  constant  of  two  orders  of  magnitude  would  be  considered  large, 

while  n  and  p  can  vary  over  ten  orders  of  magnitude  in  the  same  distance. 

2 

Thus  the  term  involving  VD  will  be  much  less  than  the  V  term.  Conse¬ 
quently,  the  Vn  •  VD^  and  Vp  •  VD^  can  usually  be  ignored  in  the  cal¬ 


culation. 


Similarly,  the  spatial  variation  of  x^  and 
lifetimes  of  holes  and  electrons  is  in  general 


T  will  be  ignored. 

.  ,  (72) 

given  by, 


The 


1_ 

T 


(16) 


where  x^ ,  x^.  •••  are  the  various  individual  radiative  and  non-radia- 

tive  lifetimes.  If  r  is  largely  determined  by  the  lifetimes  associated 
with  the  implant,  then  it  will  be  a  function  of  depth.  The  quantum 

efficiency  and  radiative  lifetime  in  p-type  GaAs  have  been  studied  by 

(73) 

Vilms  and  Spicer  .  At  77°K  they  found  bulk  quantum  efficiencies  from 

5  to  28%.  This  means  that  5  to  28%  of  the  recombination  radiation  went 

into  radiative  emission  in  the  band  which  occurs  from  0.1  ev  below  the 

band  edge  to  the  band  edge.  While  they  did  not  make  measurements  below 

77°K,  Biard,  in  a  private  communication  to  Vilms  and  Spicer,  estimated 

the  efficiency  at  30%,  substantially  independent  of  temperature.  It  is 

not  clear  whether  Biard  intended  that  to  mean  down  to  4.2  or  not.  For 

lower  implant  doses  and  higher  anneal  temperatures,  the  variation  of  the 

implant  profile  density  over  the  background  Impurity  density  is  expected 

to  be  less  than  two  orders  of  magnitude.  Tills  observation  is  based  on 

13  -2 

implant  doses  of  less  than  5  x  10  cm  ,  anneal  temperatures  of  850°  or 

15  -2 

greater,  background  impurity  concentration  of  1  x  10  cm  ,  and  the 
implant  profiles  for  Mg^^.  The  estimated  change  in  the  lifetime 
associated  with  radiative  recombination  through  the  implant  species  over 
about  the  first  .7  of  a  micron  is  one  and  one-half  orders  of  magnitude. 
Using  Biard' s  30%  figure  means  that  70%  of  the  recombination  paths  are 
non-rudiative.  These  are  assumed  to  be  independent  of  the  depth  depend¬ 
ence  of  the  implant  species.  Only  about  one-half  of  the  remaining  30%  of 
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the  recombination  paths  is  observed  to  come  from  the  recombination  path 
involving  the  implanted  species.  Thus  the  lifetime  is  a  relatively  weak 
function  of  depth  and  is  assumed  to  be  a  constant  independent  of  depth. 
Having  a  constant  lifetime  greatly  reduces  the  length  of  computer  time 
required  to  solve  the  diffusion  equation. 

The  diffusion  equation  will  also  be  restricted  to  one  dimension. 

This  assumes  a  wide  parallel  beam  of  electrons  to  be  incident  on  the 
surface  of  the  crystal.  "Wide"  is  taken  to  mean  several  times  the  dif¬ 
fusion  length  of  the  electron/hole  pairs.  The  dimensions  of  the  crystal 
are  much  larger  than  the  range  of  the  incident  electrons  and  the  ambi- 
polar  diffusion  length  of  the  generated  carriers.  This  makes  the  crystal 
appear  semi-infinite.  Charge  balance  or  electrical  neutrality  is  also 
assumed.  This  follows  the  discussion  in  McKelvey.  Since  the  experiment 
is  to  be  run  at  very  low  temperatures  (4.2°K-20oK)  the  thermal  equili¬ 
brium  values  of  electron  density  and  hole  density  are  almost  zero  com¬ 
pared  to  the  excess  carrier  density.  Thus  the  excess  hole  and  electron 
densities  are  equal.  The  experiment  is  also  conducted  in  steady-state 
so  that  all  of  the  time  derivatives  are  zero.  The  steady-state  assump¬ 
tion  also  makes  g  =  g  and  x  *  t  .  Equations  (14)  and  (15)  come  down 
p  n  p  n  ^ 

to : 


D 

P 


d26n 


dx 


„  d6n  . 

y  E  — - y  on 

p  int  dx  p 


dE 


int 

dx 


+  8n 


6n 
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„  d  5n  .  „  d6n 

D  — —  +  p  E.  — r- 
n  2  n  int  dx 
dx 


dE 


+  y  5n 
n  dx 


int 


+  g 


is. .  o 

n  in 


(16B) 


in  these  equations,  (5  n  is  the  excess  electron  density  over  the  thermal 

equilibrium  value  and  E, 

int 


is  the  internal  electric  field.  Multiplying 


terms. 


(16A)  by  bn  aiuJ  (16b)  by  p  and  adding  eliminates  both  of  the  Eint 


Dividing  by  5n  (p^  +  p^)  yields, 


D  p  +  D  U  ,2 . 

. JL21 _ 51 JJ.  iL±!l  _  ill  .  „  =  n 

U  +  U  .2  in  ®n 
a  p  ax 


(17) 


Dp  +  D  U 

Defining  D*  -  — as  the  ambipolar  diffusion  constant  results 
n  p 

in  ambipolar  transport  equation. 


n*  d  5n  6n  n 

D*  - -  -  —  +  g  =0 

,  2  Til  n 
dx 


(18) 


This  equation  is  the  same  even  if  D*  and  are  functions  of  depth. 

L(V)  Curve 

Once  the  excess  carrier  density  is  found  by  solving  equation  (18), 
the  luminescent  emission  can  be  found  by  integrating  the  rate  of  radiative 
recombination  over  the  depth  of  the  crystal.  The  rate  of  radiative 
recombination  is  proportional  to  the  number  of  states  through  which  the 
excess  carriers  can  recombine  and  to  the  number  of  excess  carriers.  The 
number  of  states  available  for  recombination  for  the  free  to  bound  Mg 
transition  is  the  number  of  Mg  ions  with  holes  bound  to  them.  As  the 
excess  carriers  recombine,  photons  are  emitted.  These  are  reabsorbed  as 
the  radiation  approaches  the  surface.  Assuming  an  exponential  absorp¬ 
tion  of  the  recombination  radiation,  the  luminescence  from  the  sample  at 
one  particular  electron  beam  voltage  can  be  found  by  summing  the  number 
of  photons  that  reach  the  surface  from  those  emitted  at  various  depths 
in  the  crystal.  This  integral  takes  the  form, 

I.(V)  =  BN(x)  5n(x)  exp  [-ax]  dx  (19) 
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where  B  is  a  constant  that  depends  on  the  geometry  of  the  experiment  and 
the  transition  matrix  element  for  recombination,  N(x)  is  ihe  density  of 
available  states  for  recombination,  and  a  is  the  absorption  coefficient 
for  the  resulting  radiation.  In  the  evaluation  of  L(V),  0  is  arbitrarily 
set  at  5  microns.  At  this  point  the  excess  carrier  density  is  down 
several  orders  of  magnitude  from  its  peak  value.  Although  N(x)  is  the 
density  of  implant  ions  available  for  excess  carriers  through  which  to 
recombine,  it  is  assumed  that  this  density  is  proportional  to  the  implant 
profile,  any  multiplying  constant  being  absorbed  into  B. 

Calculating  the  implant  profile  from  the  L(V)  curve  is  the  basic 
motivation  for  this  line  of  research.  For  this  effort,  though,  the 
Implant  profiles  are  going  to  be  assumed  and  the  calculated  L(V)  curves 
compared  to  the  measured  ones.  Thus  an  implant  profile  must  be  measured 
by  some  other  technique  or  calculated.  The  most  common  method  of  calcula¬ 
tion  Is  to  use  LSS  theory^* ^  to  determine  an  implanted  impurity  profile. 

As  the  implanted  ions  penetrate  Into  the  target,  they  lose  energy 
by  two  principal  mechanisms.  The  first  energy  loss  mechanism  is  due  to 
elastic  collisions  of  the  ion  with  the  nuclei  of  the  substrate.  These 
collisions  account  for  the  angular  scattering  of  the  ions  with  very 
'.ttle  energy  loss.  The  second  energy  loss  mechanism  is  an  inelastic 
interaction  between  the  ion  and  the  electrons  of  the  substrate  atoms. 

This  mechanism  accounts  for  the  energy  loss  of  the  ions  with  very  little 
angular  deflection.  These  energy  loss  mechanisms  are  described  by 
differential  cross  sections.  The  differential  scattering  cross  sections 
are  determined  by  the  potential  between  the  ion  and  thr  substrate  and  by 
the  nature  of  the  interaction.  LSS  theory  uses  the  Thomas-Fermi  potential 


based  upon  a  classical  statistical  model  for  the  atom.  The  LSS  results 


are  given  in  terms  of  the  mean  and  various  moments  of  a  Gaussian  distri¬ 
bution  function.  The  projected  range,  standard  deviation,  and  moment 
ratios  are  available  in  tables  tor  various  ions,  energies,  and  sub¬ 
strates.  The  available  measurement  techniques  have  been  discussed  in 
the  overview  section  of  chapter  I. 

Computer  I’.valuat  ion 

Two  programs  have  been  written  for  use  in  solving  the  diffusion 
equation.  One  actually  solves  the  equation  while  the  other  is  used  to 
evaluate  that  solution  found  by  analytically  solving  the  diffusion 
equut ion  (18) . 

The  program  that  is  used  to  solve  the  diffusion  equation  is  called 
LVCRV  for  l.(V)  curve  versus  voltage.  Given  the  depth-dose  curve  from 
EPIGA  in  the  form  of  a  data  deck,  LVCKV  fits  the  depth-dose  curve  with  a 
least  squares  approximation  by  cubic  splines  with  variable  knots.  it 
prints  out  the  voltage  being  run,  the  depth-dose  curve  from  EPIGA  and 
the  least  squares  depth-dose  curve  evaluated  at  the  same  points  as  in 
EPIGA.  The  impurity  profiles  are  computed  in  a  subroutine  called  LSS. 
There  are  provisions  for  handling  five  different  profiles  of  any  form. 
Frequently,  LSS  profiles^^  (offset  Gaussians)  are  used.  The  diffusion 
equation  Is  then  solved  for  the  excess  carrier  density.  There  Is  a 
problem  in  that  the  boundary  conditions  ate  not  specified  at  the  begin¬ 
ning  of  the  problem.  Thus  one  of  the  boundary  conditions  must  be  arbi¬ 
trarily  chosen.  The  equation  Is  then  solved.  Aa  the  solution  progresses, 
the  excess  carrier  density  Is  checked  tor  ’reasonableness."  By  "reason¬ 
ableness"  Is  meant  that  the  solution  does  not  go  negative,  nor  increase 

18 


without  bound.  The  solution  should  also  transition  into  a  decaying 
exponential  function  after  the  depth-dose  curve  goes  to  zero.  11  the 
solution  does  not  meet  these  requirements,  the  initial  condition  is 
corrected  and  the  solution  tried  again.  if  the  solution  is  not  found  in 
titty  tries,  that  ini  orm.it  ion  is  printed  out  along  wiLh  the  last  attempt. 
The  next  depth-dose  curve  is  read  in  and  the  process  repeated  until  ail 
ot  the  depth-dose  curves  have  been  run.  Finally,  the  L(V)  curve  is 
calculated  using  the  excess  carrier  densities  and  'he  impurity  profiles 
previously  calculated.  Each  point  on  the  !,(V)  curve  is  calculated  from 
equation  (19)  expressed  in  the  form  ot 

Sid 

l.(V  -  ^  LSSD(k)  *  LXN(k)  *  EXT  | -ALPHA  *  X(k)  ]  .  (20) 

k  - 1 

LSSD(k)  fs  the  impurity  density  at  depth  X(k),  EXfi(k)  is  the  excess 
carrier  density  at  the  same  depth,  Al.l’HA  Is  an  absorption  coefficient 
and  is  the  value  for  k  that  corresponds  to  X(K^)  =  V. 

In  order  to  remove  variations  in  the  experimental  data  from  point 
to  point  the  experimental  curves  are  normalized  to  an  assumed  uniformly 
distributed  impurity.  When  the  analytical  curves  are  normalized  equa¬ 
tion  (20)  is  evaluated  for  a  uniformly  distributed  impurity.  Each  point 
on  the  normalized  L(V)  curve  is  computed  by  dividing  the  result  of 
equation  (20)  for  the  implanted  impurity  by  the  result  for  a  uniformly 
distributed  impurity.  After  all  of  the  points  have  been  calculated,  the 
I.(V)  points  are  punched  on  cards  to  he  plotted  by  another  program, 

I’l.OTI.V. 

The  diffusion  equat  ion  Is  solved  using  the  Gear^^  ordinary  dilfer- 
enlial  equation  solver  package.  The  LVCKV  program  is  capable  of  handling 
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both  and  as  functions  of  depth.  This  is  the  primary  reason  for 
using  this  approach.  Since  it  takes  LVCRV  a  very  long  time  to  execute 
even  for  constant  and  p  (more  than  five  (5)  times  program  CHECK), 
another  approach  was  developed  that  would  execute  faster. 

Tliis  second  approach  is  called  CHECK.  It  utilizes  the  fact  that 
the  diffusion  equation  can  be  solved  exactly.  The  solution  to  equation 
(18)  is 

L  S  L 

6n(x)  =  N(o)  [cosh(.\/L  )  +  — rsinh (X/L  )]  -  sinh(^: — -^)g  (y)dy  (21) 

n  u7'  ii  o  Ln  n 

where  L  =  diffusion  length  =  ►  ’b*!  and 
n  n 

L 

N(o)  =  initial  condition  =  (/  expl-  y/I.  ]g  (y)dy)/(l  +  L  S/D*) 

D  *  o  n  ii  n 

This  solution  can  be  verified  by  substituting  into  equation  (18). 

Program  CHECK  Is  similar  in  many  respects  to  LVCRV.  The  depth-dose 
curves  are  read  into  the  computer  in  the  same  manner,  the  curves  are  fit 
by  the  same  subprogram,  and  they  both  use  the  same  criteria  for  a  "rea¬ 
sonable"  solution.  There  are  differences.  The  first  difference  is  that 
tiie  initial  condition  can  now  be  calculated.  Unfortunately,  the  solution 
is  very  sensitive  to  the  initial  value.  Thus  the  initial  condition  may 
have  to  be  corrected  in  order  to  get  a  viable  solution.  The  required 
changes,  though,  are  very  small,  usually  in  the  tenth  to  twelfth  decimal 
place.  The  solution,  6u(x) ,  also  has  problems  when  x  gets  greater  than 

a  tew  times  I.  .  The  subtraction  takes  place  between  two  almost  Identical 
n 

numbers  that  are  very  large,  to  get  a  comparatively  small  6n(x).  Trunca¬ 
tion  error  causes  \ n ( x )  to  oscillate  badly  under  these  circumstances. 

In  order  to  avoid  tills  condition,  the  decaying  exponential  solution  is 
substituted  for  6n(x).  Thin  solution  matches  equation  (21)  at  the  depth 
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at  which  the  depth-dose  curve  lias  its  last  value  prior  to  going  to  zero. 
The  solution  is  checked  to  make  sure  it  has  peaked  and  started  down 
before  the  exponential  solution  is  forced  on  iL. 

The  <Sn(x)  values  calculated  in  CHECK  are  printed  and  stored  in  per¬ 
manent  files.  An  additional  program,  LVCRVA,  accesses  the  permanent 
file,  calculates  the  L(V)  points  as  is  done  in  LVCRV,  scales  the  L(V) 
curves  for  plotting  and  plots  them  on  the  Cal-Comp  plotter  or  other 
plotting  device. 

Results 

As  a  confirmation  of  the  correctness  of  the  results  of  CHECK  and 
LVCRV,  L (V)  curves  were  calculated.  Figure  (10)  shows  the  results.  The 
top  curves  are  for  LVCRV  and  the  bottom  curves  from  CHECK/LVCVRA.  These 
curves  are  for  various  Implant  profiles.  The  bottom  curve  in  each  case 
is  for  a  uniformly  distributed  profile.  The  next  two  curves  are  L(V) 
curves  for  two  types  of  L.SS  profiles.  The  top  two  curves  are  the  middle 
two  curves  normalized  to  the  uniformly  distributed  profile  as  described 
earlier.  The  L(V)  curves  are  the  same  in  each  case.  In  fact,  both 
programs  give  the  same  results  for  6n(x)  to  the  five  significant  decimal 
places  printed  out. 

Since  neither  D*  t  or  S  are  known  in  advance,  a  parametric  study 
n 

was  performed  to  see  how  L(V)  changes.  A  scaling  relationship  exists 
between  the  parameters  that  substantially  reduces  the  number  of  runs 
required.  Consider  the  diffusion  equation  (18)  again. 
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—  +  gH(x)  -  0 


1)  ^ 
n  dX  | 


X  =  0 


Let  y  be  the  new  solution. 


and  let  S  be  expressed  in  units  of  /d*/t  . 

n 

The  diffusion  equation  becomes 


AD*^*  -  Ax  +  g  (x)  =  0  AU 

.  2  T  °n  n  dx 

dx  n 


x-0 


=  S/A/ (1/A)  y  x=Q  =  ASy 


x=0 


(22) 


Let  5n  =  Ay.  Then  (22)  becomes 

6n 


D*. 


d2fin 


dx* 


in 


+  gn  ( x)  =  0 


n  diSn  . 

Dnj —  =  Sfin 

dx  x=0 


x=0 


(23) 


This  is  the  same  equation  as  before.  If  I)*  is  multiplied  by  A,  di¬ 
vided  by  A  and  S  unchanged,  then  the  solution  6n  is  divided  by  A.  Since 
the  I.(V)  curves  are  arbitrary  within  a  multiplicative  constant,  they 
appear  unchanged.  The  ptiysical  significance  is  that  the  L(V)  curves 

scale  on  the  diffusion  length  given  by  L  •  /b*T  and  not  on  D*  or  6 

n  n  n 

individually.  This  also  means  that  without  an  independent  measurement 

of  ))*  or  i  ,  fitLing  the  L(V)  curve  does  not  indicate  a  value  of  either. 

This  scaling  was  verified  by  running  two  programs  with  different  D*  and 

i  ,  but  witli  the  same  L  and  S. 
n  n 

The  parametric  study  of  the  L(V)  curves  consists  of  runs  where  L  , 
S,  ALl’llA  (the  absorption  coefficient  for  the  emitted  luminescence)  and 
the  impurity  profile  are  varied.  The  L(V)  curves  are  the  result  of 
dividing  the  L(V)  curve  for  the  implant  by  the  L(V)  curve  for  a  uni¬ 
formly  doped  impurity.  This  is  done  on  a  point  by  point  basis  at  each 
voltage  where  the  L(V)  curve  is  evaluated.  This  method  of  presenting 
the  L(V)  curves  is  dictated  by  experimental  considerations  where  it  is 
necessary  to  minimize  variations  from  run  to  run  by  calibrating  the 


(a)  Diffusion  Length,  I. 


8 


(b)  Normalized  Surface  Recombination 
Velocity,  S 
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system  to  a  constant  line.  In  figure  (11)  a  series  of  runs  with  various 
combinations  of  L^,  S,  ALPHA,  and  impurity  profiles  is  shown.  The  peak 
values  of  all  these  curves  are  adjusted  to  400  so  that  conclusions  based 
on  relative  amplitudes  between  curves  are  meaningless. 

The  upper  lett  hand  figure  shows  how  L(V)  varies  as  is  changed 
over  a  factor  of  30.  The  implanted  impurity  profile  is  the  profile  pre¬ 
dicted  by  LSS  theory for  120  keV  Mg  ions  implanted  into  GaAs . 

As  expected,  the  curve  becomes  more  smeared  out  as  L  increases.  This 

n 

reflects  a  loss  of  detail  in  the  implant  profile.  It  also  demonstrates 

what  was  obvious  from  the  beginning:  the  shorter  the  diffusion  length, 

the  more  detail  of  the  impurity  profile  is  indicated  in  the  L(V)  curve. 

The  upper  right  figure  presents  L(V)  curves  with  S  as  a  parameter. 

As  the  surface  recombination  velocity  increases  there  is  more  non- 

radiative  recombination  at  the  surface  so  the  initial  parts  of  each  of 

the  excess  carrier  density  curves  (one  for  each  beam  energy)  are  reduced. 

This  shifts  the  peak  in  the  excess  carrier  density  deeper  into  the  GaAs. 

Now  lower  voltages  are  probing  areas  where  higher  voltages  were  probing 

before.  The  peak  of  the  L(V)  curve  shifts  to  the  left  as  S  increases. 

The  L(V)  curve  falls  faster  In  the  end  because  that  is  what  happens  to 

the  excess  carrier  density  as  S  increases.  The  curves  are  computed  for 

L  =  .104  micron, 
n 

The  figure  in  the  lower  left  hand  corner  shows  the  variation  of 
L(V)  with  ALPHA  as  a  parameter.  A  higher  ALPHA  reduces  the  luminescence 
more  from  deeper  parts  of  the  GaAs  to  a  greaLer  degree  than  that  from 
the  surface.  This  effectively  reduces  the  probing  effect  of  the  electron 
beam.  Hence,  higher  voltages  are  required  to  probe  deeper  into  the  GaAs 


to  find  the  implant  peak.  The  slope  of  the  initial  parts  ol  the  curve 


are  the  same;  this  indicates  there  has  not  been  much  change  in  this 
area.  The  overall  reduction  of  the  leading  edge  is  due  to  an  increase 
in  the  peak  of  the  L(V)  curve  and  the  subsequent  overall  reduction  that 
was  done  to  scale  the  plot.  The  fact  that  the  peak  and  tail  of  the  L(V) 
curve  is  larger  for  higher  ALPHA  is  at  first  surprising,  since  deeper 
effects  are  attenuated  more.  The  reason  is  that  tlie  L(V)  curve  is 
actually  the  implant  line  divided  by  a  uniformly  doped  impurity.  The 
luminescence  from  the  uniformly  doped  impurity  is  reduced  more  by  the 
higher  ALPHA  than  the  implanted  species.  Thus  the  L(V)  curves  generated 
by  the  ratio  of  implant  curve  to  doped  curve  stay  higher  for  the  higher 
ALPHA. 

The  last  set  of  curves  in  figure  (11)  illustrates  the  variation  of 
the  L(V)  curves  with  Implantation  profiles.  The  impurity  profiles  are 
LSS  profiles  for  Mg  ions  implanted  at  bO,  90  and  120  keV  into  GaAs.  The 
diffusion  length  is  0.104  microns.  As  expected,  the  L(V)  curves  for 
deeper  implants  peak  at  higher  electron  beam  energies.  The  general 
shape  of  the  curves  are  about  the  same. 

In  the  beginning  of  tills  section  it  was  shown  that  variations  in 
the  diffusion  constant  and  recombination  time  do  not  change  the  shape  of 
tiie  L(V)  curve  as  long  as  the  diffusion  length  remains  constant.  Next, 
several  paraim  trie  studies  were  made  to  see  how  the  L(V)  curves  changed 
as  tiie  diffusion  length,  surface  recombination  velocity  and  absorption 
coefficient  were  changed.  If  these  changes  are  all  taken  together,  it 
is  evident  that  some  faiiiy  complex  changes  can  occur  in  the  L(V)  curves. 
Figure  (lid)  shows  tiie  possibility  of  distinguishing  between  different 
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profiles.  This  suggests  that  cathodoiuininescenee  might  be  useful  as  a 
profiling  technique  on  either  an  absolute  basis  where  the  profile  is 
measured  directly  or  on  a  qualitative  basis  where  two  samples  are  com¬ 
pared  to  see  if  there  are  any  differences  between  them. 

The  next  section,  which  is  the  last  section  of  the  chapter,  summa¬ 
rizes  the  calculations  necessary  to  compute  the  L(V)  curve.  That  section 
is  provided  as  a  brief  review  of  the  chapter  and  as  a  quick  reference 
for  the  calculation  procedure. 

Summary  of  Calculation 

Figure  (12)  is  a  conceptual  diagram  of  how  the  calculation  of  the 
luminescence  versus  beam  voltage  (L(V))  proceeds.  As  indicated  in  the 
introduction  to  this  chapter,  the  calculation  consists  of  three  main 
parts.  The  first  part  is  the  Monte  Carlo  calculation.  This  is  charac¬ 
terized  by  module  (2).  In  this  block  the  depth-dose  curve,  which  is  the 
rate  of  energy  loss  from  the  electron  beam  with  respect  to  depth,  is 
calculated.  The  calculation  is  done  by  a  computer  program  called  EPIGA. 
Two  arbitrary  parameters,  n,  the  screening  parameter  for  the  nuclear 

charge  by  the  orbital  electrons  and  AS  ,  the  initial  step  size  for  the 

o 

Monte  Carlo  simulation,  are  adjusted  to  give  the  best  agreement  with 
Wittry  and  Kyser's  data  for  the  backscattering  fraction  of  a  30  keV 
electron  beam  penetrating  into  GuAs.  This  gives  good  agreement  with 
Martinelll  and  Wang's  data  for  electron  beam  penetration  through  thin 
films  of  GaAs. 

Next,  the  depth-dose  curves  are  used  as  inputs  to  block  (3)  which 
calculates  the  excess  carrier  density.  This  is  the  result  of  solving 
the  diffusion  equation.  In  this  calculation  both  the  diffusion  constant 
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and  recombination  time  are  lieLd  constant.  For  depths  beJow  tile  depth  at 
which  Lite  depth-dose  curves  to  to  zero,  the  excess  carrier  solution  is 
forced  to  a  decaying  exponential.  The  implementation  of  Lhis  solution 
is  done  in  two  ways  depending  on  which  of  Lhe  two  methods  is  used  to 
calculate  t  lie  excess  carrier  density. 

file  first  method  solves  the  diffusion  equation  numerically  using 
tiie  Gear  differential  equation  solver  package.  This  method  is  capable 
of  handling  both  a  depth  dependent  diffusion  constant  and  recombination 
time.  These  features  are  not  presently  used  as  the  value  of  the  calcula¬ 
tion  cannot  justify  its  cost  at  the  present  time.  In  the  second  method, 
tiie  diffusion  equaLion  is  solved  exactly  and  the  solution  evaluated 
numerically.  This  method  is  faster  than  method  one,  but  cannot  handle  a 
depth  dependent  diffusion  constant  or  recombination  time.  As  explained 
earlier  the  program  to  implement  the  first  method  is  called  LVCRV  and 
the  second  Is  called  CHECK.  The  solutions  for  either  method  agree  to 
within  the  arcu.'acy  of  the  cor  outer. 


Fig.  12  Conceptual  Diugram  of  How  the  Luni  lenccncc  Calculation  Proceeds 
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In  module  (4),  the  L(V)  curve  la  calculated  following  equation 
(20).  This  requires  the  excess  carrier  density  and  the  implantation 
profile  be  multiplied  together  at  each  depth  into  the  GaAs .  An  expo¬ 
nential  absorption  factor  Is  assumed  for  the  recombination  radiation. 
This  module  is  contained  in  program  LVCRVA  if  program  CHECK  is  used  to 
calculate  the  excess  carrier  density.  It  is  a  part  of  LVCRV  if  that 
program  is  used. 

In  module  (5)  the  data  output  is  produced.  This  consists  of  plots 
of  the  L(V)  curves.  The  plots  are  all  scaled  so  that  they  fit  onto  the 
same  axis.  This  module  is  contained  in  LVCRVA  if  CHECK  is  used  or  in 
PI.OTLV  if  LVCRV,  is  used. 

1'he  modularization  of  all  the  separate  parts  of  the  calculation  has 
made  it  easy  to  shuffle  subprograms  back,  and  forth.  Thus  modules  (4) 
and  (5)  were  combined  to  produce  LVCRVA  and  modules  (3)  and  (4)  were 


combined  for  LVCRV.  The  isolation  of  one  module  from  another  made  this 


Ill .  The  Experiment 

This  dissertation  has  as  its  primary  purpose  the  attainment  of  some 
quantitative  agreement  between  the  variation  of  the  measured  luminescent 
intensity  as  a  function  of  electron  beam  energy  and  the  same  variation 
determined  from  a  model  of  ion  implanted  GaAs.  In  order  to  do  this,  a 
facility  needed  to  be  designed  that  could  measure  luminescence  from  a 
particular  set  of  lines.  Certain  features  needed  to  be  incorporated  in 
this  system  in  order  to  get  any  kind  of  consistent  results.  These  in¬ 
clude:  an  electron  beam  which  can  be  easily  and  repeatably  controlled,  a 
beam  current  density  which  must  be  uniform  across  the  sample,  and  an 
optical  system  which  should  be  insensitive  to  small  variations  in  posi¬ 
tion.  These  requirements  were  met  to  some  extent  by  the  system  used  by 
Walter^.  The  system  was  further  modified  to  meet  the  particular  needs 
of  this  experiment. 

Epitaxial  CaAs  was  chosen  as  the  material  into  which  to  do  the  ion 

implantation.  The  epi-layer  CaAs  should  have  a  simpler  spectrum  than 

the  substrate  material  previously  studied.  Magnesium  was  selected  as 

the  implant  species  since  it  penetrates  deeper  into  the  GaAs  than  the 

(17  20) 

common  p-type  dopants  Z.n  and  Cd  “  This  was  thought  to  give  a 

better  opportunity  for  profiling  with  cathodoluminescence.  Also,  Mg  had 
not  previously  been  studied  in  GaAs  epi-layers,  although  it  had  been 
implanted  Into  substrate  quality  GaAs. 

In  the  first  sect  ion  of  this  chapter  a  description  of  the  system 
that  was  set  up  to  meet  the  particular  needs  of  this  experiment  is 
given.  Also  Included  is  a  brief  description  of  the  operating  procedures 
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used.  The  next  section  provides  some  background  on  Mg  doped  and  Im¬ 
planted  CaAs.  The  final  section  contains  a  discussion  of  the  samples 
and  how  they  were  processed. 

Experimental  Arrangement 

The  cathodoluminescence  system  can  be  broken  down  into  four  parts. 
These  are  the  electro-optical  column,  the  vacuum  system,  the  sample 
mounting  and  cooling  system  and  the  signal  processing  system.  A  de¬ 
scription  of  each  part  of  the  system  and  how  it  interfaces  with  the 
other  parts  is  given  in  the  following  sections.  A  short  description  of 
the  operating  procedures  used  with  this  system  concludes  this  section. 

Electro-optical  Column 

A  sketch  of  the  complete  system  is  shown  in  figure  (13) .  The  elec¬ 
tro-optical  column  consists  of  the  vacuum  chamber,  the  steering  and 
focus  coils,  and  the  electron  gun  system.  The  electron  gun  is  a  B-91- 
16A35  unit  manufactured  by  Hughes  Aircraft  Company.  The  electron  gun  is 
attached  to  a  circular  eight-pin  connector  which  plugs  into  the  eight 
pins  of  the  glass  envelope  surrounding  the  gun.  A  centering  clip  is 
attached  to  the  anode  to  center  the  gun  in  the  glass  envelope  and  also 
provide  a  ground  connection  for  the  anode.  There  are  two  grids  on  the 
electron  gun.  The  grid  voltages  are  variable  but  usually  G1  is  about 
-15  volts  and  G2  is  +500  volts.  The  electron  gun  is  operated  at  pres¬ 
sures  ranging  from  5  x  10  ^  torr  to  1  x  10  ^  torr. 

After  the  electron  beam  leaves  the  gun,  it  passes  through  an  aperture 
positioned  in  the  center  of  a  45°  mirror.  Interchangeable  stainless 
steel  apertures  are  provided  with  diameters  of  5,  8,  10,  12  and  15  mil. 
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These  apertures  were  not  used  since  they  reduce  the  amount  of  current 
the  gun  can  deliver  Lo  the  target.  The  mirror  is  positioned  so  that  it 
is  possible  to  look  straight  up  the  column  to  the  samples.  The  electron 
gun  housing  is  in  a  sLiding  mount  below  the  mirror  block.  This  mount 
may  be  moved  transversely  on  a  sliding  O-ring  seal.  This  corrects  for 
misalignment  of  the  electron  gun. 

There  are  five  sets  of  magnetic  coils  that  the  beam  passes  through 
before  it  reaches  the  target.  They  are  the  alignment  coils,  centering 
coils,  focus  coil,  stigmator  coils,  and  deflector  coils.  The  coils  are 
all  connected  to  D.C.  power  supplies  capable  of  providing  at  least  200 
milliamps  into  the  coils.  All  but  the  focus  coil  have  coils  for  both 
horizontal  and  vertical  beam  adjustments.  The  alignment  coils  are  on 
the  electron  gun  housing  and  slide  with  it.  Coupled  with  manually 
sliding  the  gun,  the  alignment  coils  are  used  to  get  the  maximum  current 
and  most  uniform  spot  through  the  aperture.  This  Is  determined  by 
viewing  the  electron  beam  on  a  movable  flap  that  has  a  zinc-oxide  crystal 
mounted  on  it.  The  flap  Is  located  at  a  position  immediately  after  the 
beam  passes  through  the  aperture.  The  centering  coils  are  used  to  posi¬ 
tion  the  beam  on  the  axis  of  the  focus  coil  so  as  to  get  the  most  uniform 
spot  on  the  target.  The  current  through  the  focus  coil  is  adjusted  so 
that  the  beam  is  defoeused.  The  spot  size  is  enlarged  until  the  spot 
covers  the  entire  sample  exposed  to  the  beam.  This  is  necessary  so  that 
the  current  density  at  the  surlace  of  the  GaAs  can  be  approximately 
controlled.  Later,  It  will  be  shown  that  failure  to  control  the  current 
density  makes  a  considerable  dillereiice  in  the  spectrum.  The  stigmator 


coils  correct  for  any  stigmatism  in  the  beam.  It  was  not  necessary  to 


use  Liu*  sL  iguutor  coils.  The  de  1  1  ec  t  ion  I'oils  are  used  Lo  steer  the 
beam  onto  the  target  and  into  the  Faraday  cup.  A  switch  is  used  to 
se 1 ec  L  which  of  two  Sets  of  two  power  supplies  is  connected  to  the 
coils.  this  allows  the  beam  to  be  conveniently  switched  from  tile  sample 
to  the  Faraday  cup  and  back.  The  focus  coil  and  the  stigmator  coils  are 
located  on  a  uiicroposit  ioner  that  features  piLcii  and  yaw  tilt  adjustment 
and  two  axis  translation.  The  deflection  yoke  is  mounted  on  a  similar 
microposi t ioner  that  also  allows  for  rotation  of  the  coils. 

The  beam  current  is  measured  by  means  of  a  Faraday  cup  and  a  Keith- 
ley  414A  p icoainmeter .  The  beam  enters  the  cup  through  an  1/8  inch  dia¬ 
meter  hole.  This  is  the  same  size  hole  that  is  in  the  shield  that  holds 
the  samples  in  place.  By  detocusing  the  beam  so  that  It  covers  the 
entire  bole  and  by  assuming  the  beam  to  be  uniform  across  the  spot,  the 
current  density  can  be  adjusted  to  any  desired  value.  Any  variation  in 
the  uniformity  of  the  beam  cross-section  is  accounted  for  by  adjusting 
the  position  of  the  beam  on  the  cup  to  maximize  the  current  reading  and 
on  the  sample  to  maximize  the  luminescence  Lo  the  photomultiplier.  This 
adjustment  allowed  the  amplitude  measurement  of  the  luminescence  to  be 
repeated  to  within  10%  for  subsequent  measurements  that  were  immediately 
repeated.  The  shield  over  the  entryway  to  the  Faraday  cup  is  kept  at  the 
negative  grid  voltage  of  the  suppression  screen. 

Vacuum  System 

Figure  (13)  includes  the  major  components  of  the  vacuum  system. 

Kef or once  to  this  diagram  shows  how  the  system  components  interrelate. 

The  lw)  diffusion  pumps  allow  the  samples  to  be  changed  without  bringing 
the  gun  up  to  atmospheric  pressure. 


The  sample  mounting  Ls  the  same  as  has  been  used  in  the  past.  Four 
CaAs  samples  are  mounted  on  the  face  of  a  copper  block  that  has  been  cut 
at  a  450  angle  to  the  electron  beam  axis  and  the  horizontal  axis  of  the 
spectrometer.  The  samples  are  held  in  place  with  copper  brackets  that 
are  attached  with  screws  to  Lhe  copper  block.  The  brackets  have  an  1/8 
Inch  diameter  hole  in  them  to  allow  the  electron  beam  to  strike  the 
sample  and  the  luminescence  to  be  collected.  The  brackets  are  coated 
with  aquadag  to  minimize  reflection  from  the  electron  gun  filament.  A 
Faraday  cup  is  mounted  at  the  top  of  the  face  of  the  copper  block. 

An  Air  Products  lleli-Tran  system  is  used  to  cool  the  samples. 

Liquid  helium  flows  through  the  Lip  of  the  transfer  tube,  cooling  the 
copper  block.  The  temperature  of  Lhe  copper  block  is  monitored  with  a 
chromel  vs.  gold  .07  atomic  Z  iron  t henuoeoupl e  using  a  liquid  nitrogen 
reference.  A  Data  Precision  1500  digital  voltmeter  is  used  to  record  the 
temperature,  in  microvolts. 

S i gi  1  al  Processing 

The  luminescence  from  the  sample  is  collected  by  a  three  lens  system 
as  shown  in  figure  (13).  lhe  lens  closest  to  the  samples  has  a  focal 
Kngth  of  10  cm  and  is  positioned  approximately  this  distance  from  the 
samples.  The  other  lenses  are  In  Lhe  same  mount  with  a  fixed  separation 
of  8  cm.  They  were  positioned  near  the  spectrometer  slit  and  adjusted  to 
maximize  the  output  signal.  The  lens  closest  to  the  slit  has  a  focal 
length  ol  6.5  cm  and  the  other  lens  a  local  length  of  25  cm.  Between  the 
second  and  third  lenses  is  a  number  650  Spectracoat  Varipass  filter.  It 
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is  used  to  prevent  any  of  the  Incident  radLation  below  6500  A  from  entering 
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the  spectrometer.  Stray  light  from  the  filament  is  also  reduced  by  sur¬ 
rounding  the  electron  gun  with  a  shroud  made  from  a  black  cloth. 

The  luminescence  is  focused  onto  the  input  slits  of  a  Spex  model 
1702  Czerny-Turner  spectrometer.  A  200  micron  slit  opening  is  used.  A 
600  g rooves/mill iineter  Bausch  and  Lomb  grating  blazed  at  1.6  microns  is 
used  in  second  order  to  disperse  the  luminescence.  The  efficiency  of  the 
grating  is  estimated  to  vary  by  about  10%  over  the  range  of  wavelengths 
of  interest  here,  being  more  efficient  at  the  longer  wavelengths.  The 
luminescence  is  detected  by  an  RCA  C7000A  photomultiplier  tube  with  an  S- 
1  response.  Hie  response  changes  by  about  10%  over  the  same  wavelength 
range,  decreasing  for  longer  wavelengths.  Thus  the  two  responses  tend  to 
cancel  each  other,  making  the  response  of  the  system  relatively  flat  over 
the  frequency  range  used  here.  The  tube  fs  cooled  to  -50°C  using  liquid 
nitrogen.  The  temperature  is  maintained  at  this  value  by  a  Products  for 
Research  temperature  controller.  The  output  ol  the  detector  is  amplified 
by  a  PAR  1121  discriminator/araplifier.  While  the  discriminator/amplifier 
provides  many  other  functions  such  as  puLse  height  analysis,  window 
scaling,  and  pulse  pile-up  correction,  the  unit  is  operated  in  the  single 
mode.  in  this  mode  only  pulses  greater  than  a  pre-set  level  are  detected 
and  sent  on  to  the  next  stage.  file  pre-set  level  is  determined  by  a 
trial  and  error  procedure  that  maximizes  the  number  of  signal  counts 
while  minimizing  the  number  of  noise  counts.  Unce  this  level  Is  selected, 
it  is  not  changed  during  the  experiment. 

The  next  step  in  processing  the  signal  is  to  shape  the  pulses  so 
that  they  are  compatible  with  the  inpuL  to  the  multichannel  analyzer. 

The  Hewlett-Packard  5400A  multichannel  analyzer  used  in  this  experiment 
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requires  pulses  of  at  least  4  volLs  height  and  25  nanoseconds  widLh.  A 
Hewlett-Packard  8010A  pulser  operating  in  the  negative,  external  trig¬ 
gered  mode  Is  used  lor  pulse  shaping  and  interfacing  between  the  dis¬ 
criminator/amplifier  and  the  multichannel  analyzer. 

The  multichannel  analyzer  is  operated  in  the  MCS  mode.  In  this  mode 
the  pulses  from  the  pulse  generator  are  counted  for  a  specified  length  of 
time  either  with  an  internal  or  external  time  base.  The  number  is  stored 
in  one  channel  of  the  unit's  memory.  This  process  is  repeated  over  and 

over  until  all  of  the  channels  have  numbers  stored  in  them.  The  analyzer 

has  1024  channels  of  memory.  For  this  experiment,  the  multichannel 
analyzer  Is  operated  under  the  external  controls  of  an  interface  control 
box  (see  figure  13).  This  device  determines  the  length  of  time  in  which 
pulses  are  counted  for  each  channel.  This  controller  is  also  used  to 
drive  the  spectrometer,  thus  synchronizing  the  spectrometer  and  the 
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multichannel  analyzer.  The  controller  has  been  described  elsewhere 
and  lienee  will  not  be  repeated  here. 

The  data  stored  In  the  multichannel  analyzer  Is  read  out  on  to  paper 
tape.  The  paper  tape  is  taken  to  a  computer  terminal  site  where  the  tape 

is  converted  to  punched  cards.  The  cards  are  used  in  a  computer  program 

that  produces  an  output  via  a  Cal-Comp  or  similar  plotter.  The  plots  are 
of  arbitrary  intensity  versus  wavelength.  (It<-  abscissa  of  the  plots  is 
calibrated  using  a  mercury  lamp  to  ptovide  a  calibration  line  over  the 
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spectrum  when  the  original  GaAs  Luminescence,  is  recorded.  The  4046.6  A 
llg  line  is  used  In  second  outer  for  the  calibration  line. 
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Operating  Procedures 

In  this  section  the  procedures  used  during  the  experiment  are  ex¬ 
plained.  The  purpose  of  the  section  is,  (a)  to  simplify  follow-on  ex¬ 
periments  Lhat  use  the  equipment  and,  (b)  to  describe  the  operating 
conditions  under  which  the  data  is  taken. 

The  first  step  is  to  choose  the  four  samples  to  be  measured  and  put 
them  in  place.  The  samples  already  in  the  chamber  (assuming  there  are 
samples  in  the  chamber)  are  isolated  from  the  vacuum  pumps  by  means  of 
the  gate  valves.  That  portion  of  the  chamber  containing  the  samples  is 
isolated  from  the  electron  gun  section  and  brought  up  to  atmospheric 
pressure  by  filling  the  chamber  with  nitrogen.  The  cold  finger  is 
removed  and  Lhe  samples  changed.  The  procedure  is  reversed  to  get  the 
samples  back  Into  the  chamber. 

The  samples  are  then  cooled  using  the  Heli-Tran  system.  The  system 
temperature  Is  let  fall  to  its  lowest  level.  The  liquid  helium  flow  is 
then  reduced  to  just  maintain  this  temperature.  The  estimated  sample 
temperature  is  10°K. 

Next,  the  electron  beam  is  directed  onto  a  sample.  The  optical 
system  is  aligned  using  the  luminescence  from  the  sample  as  a  source. 

With  the  system  aligned,  the  electron  beam  is  switched  to  the  Faraday  cup 
(see  the  Electro-optical  Column  section).  The  beam  current  is  adjusted 
to  the  desired  level  (usually  r) .  0  microamps)  with  the  beam  spot  defocussed 
to  cover  the  entire  opening  in  Lhe  Faraday  cup  shield.  The  beam  is 
steered  onto  the  cup  to  give  Lhe  maximum  reading.  The  beam  is  then 
switched  back  to  the  sample  and  lhe  current  to  the  deflection  yoke  adjusted 
to  steer  the  beam  into  a  position  of  maximum  lumi nescence  at  a  particular 
Line. 


The  data  is  now  recorded  in  an  hp-5400A  multichannel  analyzer. 
Sometimes  the  entire  band  of  luminescence  from  1.3  to  1.5  eV  is  recorded. 


At  other  times  only  certain  lines,  as  for  example,  the  Mg  complex  and  the 
copper  complex,  are  recorded.  When  a.1  l  1024  channels  of  the  analyzer  are 
filled,  the  data  is  transferred  to  punched  paper  tape  and  then  to  punched 
cards.  These  cards  can  be  used  in  a  program  for  plotting  the  data. 

Plots  of  the  luminescence  are  then  analyzed  on  a  Dupont  310  line- 
shape  simulator.  This  device  takes  up  to  six  waveshapes  and  combines 
them  together  to  get  a  composite  line.  This  composite  line  can  then  be 
compared  to  the  real  line.  Thus  competing  lines  can  be  resolved  by 
adjusting  the  individual  lines  making  up  the  manufactured  composite  line 
until  it  matches  the  complex  waveshape.  The  other  lines  can  be  switched 
off  leaving  only  the  one  line  to  study.  The  device  also  determines  the 
area  under  the  curve.  This  is  the  integrated  lineshape  and  corresponds 
to  the  total  number  of  transitions  through  that  center  per  unit  time. 

The  L(V)  curves  are  calculated  by  taking  the  ratio  of  the  area  under  the 
Mg  curve  to  the  area  under  the  Cu  curve  at  each  voltage. 

Mg/GaAs  Luminescence:  Previous  Results 

The  first  reported  case  of  Mg  being  implanted  into  GaAs  seems  to 

have  been  by  Hunsperger  et  al^^.  They  implanted  Be  and  Mg  at  45  keV 

into  n-type  (n  =  1.0  x  1016/cm3)  GaAs  substrates.  The  implantation  was 

done  at  room  temperature  with  subsequent  annealing  at  temperatures  from 
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1 00-900°C .  The  implanted  GaAs  was  capped  with  a  2000  A  thick  layer  of 
SiO^  prior  to  annealing.  They  determined  that  the  p-type  layer  was 
approximately  0.4  microns,  whereas  LSS  theory  predicted  0.12  microns. 
They  concluded  that  some  ion  channeling  or  diffusion  had  occurred  and 
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Chat  it  could  not  be  assumed  that  the  implanted  layer  had  a  Gaussian 
shape  as  predicted  by  LSS  theory. 

Yu  and  Park^^  used  photoluminescence  to  study  the  luminescent 

behavior  and  the  annealing  characteristics  in  Mg  ion  implanted  substrate 

materials.  Si-doped  and  Cr-doped  semi-insulating  substrates  were  used, 

with  the  Mg  ions  being  Implanted  directly  into  the  substrates.  Emission 
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spectra  were  obtained  at  4.2K  for  implanted  doses  of  from  5  x  10  /cm  to 
15  2 

1  x  10  /cm  and  for  annealing  temperatures  of  450  to  900°C.  The  con¬ 
clusions  were  that  the  optical  ionization  energy  of  Mg  is  28  ±  2  meV  and 
that  annealing  temperatures  from  750-900°C  are  enough  to  optically  acti¬ 
vate  implanted  Mg  ions  and  to  remove  lattice  damage.  The  luminescence 
from  the  Si-doped  CaAs  showed  three  bands.  One  band  at  1.513  eV  was 
attributed  to  impuricy-exciton  complexes.  The  middle  band  that  occurred 
at  various  places  between  1.493  to  1.489  eV  was  identified  as  free  elec¬ 
trons  recombining  with  holes  bound  to  acceptors.  The  lower  band  comes 
from  electrons  bound  to  donors  recombining  with  holes  bound  to  the  same 
acceptors  as  before.  The  position  of  this  bund  was  very  much  intensity 
dependent,  centered  around  1.464  eV.  Similar  results  were  obtained  for 
the  Cr-doped  substrate  except  that  the  lowest  band  did  not  occur  at  a 
lower  annealing  temperature  (750°C). 

In  a  follow-on  paperv  '  the  temperature  dependence  of  photolumines¬ 
cence  from  Mg  implanted  GaAs  was  studied.  The  same  samples  as  in  refer¬ 
ence  (77)  were  used.  In  addition,  an  undoped,  n-type  substrate  implanted 
with  Mg  was  studied.  Again  the  major  observed  bands  were  from  free  elec¬ 
trons  recombining  with  holes  bound  to  Mg  acceptors  and  from  donor-acceptor 
pair  recombination.  The  dose  dependence  of  the  emission  characteristics, 
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the  temperature  dependence  of  emission  bands  involving  Mg,  and  the  donor- 
acceptor  pair  band  were  all  studied.  The  depth  distribution  of  Mg  accep¬ 
tors  obtained  by  successive  etchings  indicated  that  Mg  acceptors  were 
distributed  to  about  0.65  microns  for  samples  annealed  at  900°C.  The  LSS 
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depth  for  the  100  keV  implant  is  1030  A  with  a  straggling  of  460  A.  Yu 
and  Park  conclude  that  the  disparity  is  caused  by  diffusion  of  the  Mg 
ions  during  annealing. 

(79) 

In  a  more  recent  paper,  Yu  discussed  the  excitation  dependent 
emissions  in  Mg,  Be,  Cd,  and  Zn  implanted  GaAs  with  respect  to  changes  in 
temperature  and  excitation  intensity.  He  observed  the  donor-acceptor 
pair  emission  in  the  impure-compensated  crystals  shift  to  lower  energy 
with  a  temperature  increase  in  the  range  of  4-50K  while  the  donor-accep¬ 
tor  pair  emission  in  pure  crystals  behaves  in  the  usual  manner.  The 
large  shift  of  the  emission  peaks  is  thought  to  occur  in  the  impure, 
compensated  regions. 

Zolch  et  al^^  studied  implants  of  Be,  Cd,  Mg  and  Zn  in  GaAs  and 
GaAs^^P^.  For  Mg  implants  of  doses  between  10^  and  lO^/cm^  they 
observed  a  reverse  annealing  at  temperatures  between  650°C  and  750°C.  By 
this  they  mean  the  effective  mobility  decreased  for  annealing  tempera¬ 
tures  in  this  range.  They  also  observed  that  a  higher  electrical  acti¬ 
vation  and  a  smaller  reverse  annealing  effect  can  be  obtained  if  Si^N^ 
passivation  is  used  rather  than  SiC^.  Annealing  was  done  in  a  nitrogen 
atmosphere.  The  reverse  annealing  effect  was  attributed  to  outdif fusion 
of  Mg  and  not  complex  formation  since  outdif  fusion  is  reduced  with  Si^N^,. 
Annealing  above  800°C  was  required  before  the  mobility  of  the  implanted 
layer  approaches  the  bulk  value.  They  found  the  diffusion  constant  of  Mg 
to  be  (kT  in  eV) 
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D  =  4.4  x  10  4  exp  (-1.87/kT)  cm^/sec. 

Unfortunately,  they  did  not  report  their  profile  measurements  for  Mg. 

Samples 

Mg  implanted  epitaxial  CaAs  was  chosen  for  study  to  compare  experi¬ 
mental  with  theoretical  results  for  the  luminescence  versus  voltage 
(L(V))  curves.  Epitaxial  GaAs  was  chosen  since  the  spectrum  from  epi- 
layers  should  be  simpler  to  interpret  as  there  are  less  competing  losses. 
Mg  was  picked  for  the  implant  for  two  reasons.  First,  no  previous  work 
had  been  done  with  Mg  implanted,  epitaxial  GaAs.  Second,  the  LSS  pro¬ 
jected  range  (appendix  A)  of  Mg  ions  is  from  .06  to  .12  micron  for  the 
implant  energies  available  to  us.  This  is  about  the  same  depth  as  that 
of  the  peak  of  the  depth-dose  (see  fig.  3  and  5-9).  Having  the  peak  of 
the  depth-dose  curve  shift  through  the  peak  of  the  implant  profile 
increases  the  chances  of  profiling. 

GaAs  epilayers  were  implanted  at  fluences  of  1E12,  5E12,  1E13,  5E13 
and  lE14/cm^.  The  implant  energies  used  were  60,  90  and  120  keV .  The 
samples  were  then  cupped  with  a  layer  of  Si^N^,  annealed,  and  stripped  of 
the  Si^N^  layer.  Finally,  some  were  profiled  using  Hall  measurements 
coupled  with  repeated  etchings.  These  topics  in  more  detail  in  the  fol¬ 
lowing  sections. 

Growing 

The  epi-layers  were  grown  on  a  chromium  doped  GaAs  substrate  using  a 
vapor  phase  epitaxial  growth  technique.  They  were  grown  by  G.  McCoy  of 
the  Electronics  Research  Branch  of  the  Air  Force  Avionics  Laboratory  at 


Wright-Patterson  Air  Force  Base.  An  open  flow  chemical  vapor  transport 


system  employing  the  reagents  arsenic  trichloride,  gallium,  and  hydrogen 
was  used.  The  system  is  described  in  more  detail  in  reference  (81).  The 

epilayers  were  p-type  with  a  room  temperature  carrier  concentration  of 

14  -i  2 

7.94  x  10  cc  and  an  effective  mobility  of  199  cm~/v  sec  as  determined 

by  Hall  measurements. 

Implanting 

The  epilayers  were  implanted  by  J.  Ehret  of  Systems  Research  Labora¬ 
tories  using  the  ion  implantation  system  at  the  Avionics  Laboratory. 

Before  implanting,  the  samples  were  cleaned  in  methanol.  A  hot  cathode 
ion  source  of  solid  Mg  was  used  as  a  source  of  Mg.  The  Mg  ions  were 
separated  from  the  other  ions  by  a  magnetic  separator.  The  Mg  ions  were 
then  accelerated  to  the  required  energy  for  implantation  into  the  samples. 

More  information  on  the  system  is  contained  in  reference  (84). 

Table  II  provides  a  listing  of  the  GaAs  samples  available  for  study. 

A  large  number  of  samples  was  needed  since  it  was  not  known  beforehand 
which  implant  energy  and  fluence  would  give  the  best  L(V)  curve.  It  was 

felt  that  the  mid-range  doses  would  be  best  since  too  low  a  dose  decreases 

the  amount  of  luminescence  and  too  large  a  value  leads  to  impurity  banding 
(which  is  another  complication  in  the  interpretation  of  the  data).  The 
unimplanted/annealed  and  the  unimplanted/unannealed  samples  were  used  as 
references  to  see  the  effects  of  implanting  and  annealing  and  to  aid  in 
the  identification  of  the  lines  in  the  implanted  samples. 

Anneal lng 

Table  II  shows  the  conditions  undex  which  the  samples  were  annealed. 
Si^N^  caps  were  used  since  at  the  time  those  were  giving  the  most  repeatable 


Table  II 
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results''  .  The  850°  anneal  for  fift  een  minutes  was  somewhat  arbi¬ 
trarily  set.  It  has  been  shown^^’  that  annealing  at  temperatures 

greater  than  750°C  is  sufficient  to  obtain  excellent  electrical  proper¬ 
ties  and  efficient  luminescence  from  Mg  acceptors.  The  usual  anneal  time 
has  been  fifteen  or  twenty  minutes  .  The  samples  to  be  annealed 

were  covered  with  a  pyrolytically  grown  film  of  Si^N^  approximately  1000 

O 

A  thick.  They  were  then  placed  in  a  spectrosil  boat  and  annealed  at  the 
required  temperature  in  flowing  hydrogen.  After  the  samples  cooled,  the 
caps  were  removed  using  hydrofluoric  acid  with  trichloroethylene,  acetone 
and  de-onlxed  water  rinses.  ft  was  decided  to  remove  the  caps  in  order 
for  the  peak  of  the  electron  beam  to  probe  as  deeply  as  possible  into  the. 
sample . 

1  l_e ctri cal  H r ol  i  1  lug 

-2 

After  the  L(V)  curves  were  obtained,  the  3E13  cm  at  120  keV  and 
_  > 

the  5 K 1 3  cm  “  at  90  keV  samples  were  profiled  using  a  Hall  system  with 
repeated  chemical  etchings.  The  implant  density  was  determined  about 

O 

every  L50  to  200  A.  The  volume  concentration  of  charge  carriers  and  the 

mobility  can  tlu-n  he  calculated  as  a  t unction  of  depth.  The  profiles 

were  measured  by  V.  K.  Yeo  of  the  Electronics  Research  Branch.  Figure 

(14)  shows  the  resuLts  for  the  120  keV,  5K13  cm  *"  and  the  90  keV,  5E1 3 
-2 

cm  "  Implant.  The  points  show  considerable  scatter  compared  to  other 
profiles  on  similar  samples.  'The  implications  of  this  will  be  developed 
further  in  the  chapter  on  results.  With  this  equipment  and  these  samples, 
the  required  data  can  be  acquired.  This  data  Is  reported  in  the  next 
chapter . 
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IV.  Results 


This  chapter  contains  both  the  results  obLained  using  the  system 
described  in  the  last  chapLer  and  the  comparison  of  the  experimental  and 
analytical  l.(V)  curves.  First,  the  spectrum  from  Mg  ion  implanted  GaAs 
is  analyzed.  Various  Lines  between  1.3  and  1.5  eV  are  identified.  It  is 
shown  how  depth  resolved  cat hodoiuminescence  can  be  used  to  identify 
lines  by  determining  where  the  concentration  of  the  impurity  is  maximum. 
Next,  a  study  is  conducted  of  the  spectral  variation  with  current  density. 
A  possible  explanation  for  this  variation  is  given.  The  reason  that  this 
variation  is  important  to  depth  resolved  cathodoluminescence  is  discussed 
and  the  proper  a  t  ion  taken  to  minimize  its  impact  is  given.  The  L(V) 
curves  are  presented  in  the  third  section.  In  the  final  section  the 
analytical  and  experimental  L(V)  curves  are  compared.  The  physical  para¬ 
meters  chosen  for  the  analytical  curves  are  justified  based  on  other 
workers'  results. 

Spec  L rum  ldcnl i f icat 1  on 

Figure  (15)  shows  the  spectrum  obtained  from  a  GaAs  epllayer  im¬ 
planted  at  120  keV  to  a  dose  ot  5E13  cru  with  Mg  ions.  The  spectrum  was 
obtained  with  a  15  keV  beam  of  electrons  at  a  5  microamp  Faraday  cup 
current.  (The  total  cup  current  will  usually  be  specified  instead  of  the 
current  density.  In  all  cases,  the  beam  of  electrons  was  defocused  to 
cover  the  entire  entrance  aperture  to  the  cup  so  that  when  the  total  cup 
current  Is  the  same,  so  is  the  current  density  for  a  uniform  beam  is  als* 
identical.)  The  spectrum  has  not  been  corrected  for  the  response  of  the 
detector,  although  the  response  should  be  approximately  flat  over  this 


I requeucy  range. 
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In  general,  there  are  four  or  five  bands  of  luminescence.  The  1.560 
eV  line  has  been  attributed  either  to  copper  or  to  an  As 

vacancy  bound  to  an  acceptor .  a  1.350  eV  line^^^  has  been  identi¬ 
fied  but  it  is  too  far  from  1.36  eV  to  be  considered.  Chiang  and  Pearson 
found  that  the  peaks  could  be  enhanced  with  Cu  diffusion  or  annealing  in 
Amersil  quartz,  but  were  absent  in  samples  annealed  in  Spectrosil  quartz. 
The  peaks  at  1.325  and  1.288  eV  are  longitudinal  optical  phonon  replicas 
of  the  1.360  eV  line.  References  (87)  and  (89)  show  that  the  peaks  they 
studied  follow  a  Poisson  distribution.  The  ratio  of  the  peaks  in  this 
study  approximately  follows  that  ratio  (see  table  III).  Thus  the  line  at 
1.36  eV  is  the  same  line  as  the  copper  line  identified  by  Queisser  and 
Fuller. 


fable  111 

Ratio  of  1.288  eV  and  1.324  eV  Peaks  to  1.360  eV  Peak 

Peak  Kxper imental  Ratio  Poisson  Ratio 

1.325  eV  .424  .43 

1.288  eV  .118  .09 

The  identity  of  the  line  at  1.408  eV  is  not  clear.  It  is  present  in 
the  unimplanted/capped  and  annealed  sample  and  all  of  the  implanted/ 
capped  and  annealed  samples.  It  is  absent  in  the  unimplanted/unannealed 
sample.  The  line  Is  not  associated  with  a  vacancy  complex  as  the  anneal¬ 
ing  process  would  be  expected  to  reduce  its  Intensity.  A  line  at  1.408 
eV  has  been  associated  with  Mn  or  Pb^^.  The  line  at  1.400  eV  is  con¬ 
sidered  a  transverse  acoustic  (TA)  phonon  replica  of  the  1.408  eV  line. 
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In  addition  to  its  energy  spacing.  It  is  not  present  in  the  spectrum  when 
the  1.408  eV  is  not  there.  Other  samples  measured  have  a  weak  LO  phonon 


replica  from  the  1.408  eV  line.  This  structure  is  similar  to  that  shown 
in  references  (87)  and  (90)  for  Mn.  Otherwise  there  does  not  seem  to  be 
any  particular  reason  to  choose  Mn  or  l’b.  Since  the  spectrum  of  the 
un implanted/unannealed  samples  does  not  have  the  1.408  eV  line,  it  must 
originate  in  the  implantation  or  annealing  process.  Examining  the  ratio 
of  the  L.408  eV  line's  intensity  to  the  copper  line's  intensity  as  the 
beam  voltage  (and  hence  probing  depth)  is  varied  should  give  some  idea 
whether  the  intensity  originates  primarily  near  the  surface  or  from 
deeper  into  the  epilayer.  This  would  IndicaLe  whether  the  impurity 
diffused  into  the  epilayer  from  tin-  surface  or  the  substrate.  Table  IV 
shows  the  ratios  for  two  samples.  the  conclusion  is  that  the  impurity 
diffused  into  the  epilayer  from  the  1  rout  surface  during  either  the 
capping  or  annealing. 


Table  IV 

Ratio  L.408  eV/Cu  Line 


Beam  energy  (keV  at  5  microamp)  Sample  1  Sample  2 

5  .453  x 

7.5  .247  x 

10  .178  .25 

12.5  .117  .216 

15  .104  .132 

_2 

Sample  1:  90  keV  Mg/lElt  cm 

Sample  2:  90  keV  Mg/lKli  cm 


The  complex  around  1.49  eV  is  m^Je  up  of  lines  at  1.48b,  1.490,  and 

1.492  eV.  The  1.490  eV  peak  seems  likely  to  be  due  to  Zn  since  Ashen  et 
(91 ) 

al  reported  finding  zinc  in  over  95Z  of  their  vapor  phase  samples 

with  the  other  5%  probably  having  it.  They  report  the  peak  at  1.4888  eV 
buL  pulling  of  tlie  peak  by  the  much  stronger  peak  at  1.492  eV  seems 
likely.  The  Zn  line  is  present  in  alL  of  the  samples  including  the  unan- 
nealed/unimplanted  and  the  anneu 1  ed/un  implanted .  The  1.492  eV  line  is 
assigned  Lo  Mg.  The  line  energy  to  tour  significant  digits  is  1.4918  eV. 
The  assignment  in  Ashen,  et  al  is  1.4911  eV  ( '-  .0003  eV)  .  The  uncer¬ 
tainty  in  the  position  of  the  peak  is  about  '1  channel  (±  .0004  eV) .  The 
thermal  distribution  of  electrons  in  the  conduction  band  subtracts  !jkT 
from  the  energy  of  the  line  (see  Ashen).  This  brings  the  energy  assign¬ 
ment  well  into  agreement.  This  line  is  absent  in  both  tin.*  unannealed/ 

unimplanted  and  uiuieaied/un implanted  samples. 

(91)  (91) 

The  line  at  1.486  eV  could  be  Si  ,  Cd  or  a  donor-acceptor 

pair  with  the  1.492  eV  line^^.  Since  the  line  at  1.486  eV  is  present 
when  tlie  1.492  eV  Ls  not,  the  donor -acceptor  pair  assignment  is  ruled 
out.  As  for  Si  or  Cd ,  Lhe  line  is  assigned  to  Si.  In  the  unimplanted/ 
annealed  sample  lhe  1.486  line  is  much  stronger  than  the  1.486  line  in 
unimplanted/  unannealed  sample.  This  is  attributed  to  diffusion  of  Si 
into  the  sample  during  the  capping  and  annealing.  Also  if  the  SI  is 
coming  In  from  the  cap,  tl  silicon  line  should  grow  compared  to  the 
other  lines  as  the  electron  beam  probes  shallower  depths.  This  ls 
clearly  seen  In  the  sequence  oi  runs  shown  in  figures  (17)  and  (18).  In 
these  figures  the  1.486  eV  line  clearly  increases  compared  to  the  Zn  line 
as  the  beam  energy  decreases,  lionet  probing  shallower  depths.  Figures 
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( L  7 )  and  (18)  not.  only  support  the  Si  assignment ,  Liny  also 
the  profiling  possibilities  of  eathoilo  I  umi  neseence  . 

The  line  at  1.449  is  a  TO  phonon  replica  ol  the  line  at  1.486  eV. 

The  energy  spacing  is  .037  eV  which  is  approximately  the  1.0  phonon 

/  Q  -7  \ 

energy  .  The  lines  at  1.455  and  i .453  eV  are  also  .037  eV  below  the 

lines  at  1.492  and  1.490  eV,  respectively,  and  are  also  phonon  relicas. 

In  some  samples  the  lines  around  1.45  eV  showed  the  exact  shape  of  the 
lines  around  1.49  eV  (see  figure  16).  When  the  1.492  eV  line  is  absent 
from  the  spectrum,  so  is  the  line  at  1.455  eV. 

The  line  at  1.504  eV  is  unknown  at  present.  It  is  absent  from  the 
samples  that  were  not  implanted  with  Mg,  leading  to  the  conclusion  that 
it  is  associated  with  Mg.  On  the  other  hand,  it  is  also  missing  from 
some  of  the  Mg  implanted  samples.  It  is  9  meV  from  the  1.513  eV  line, 
which  suggests  it  might  be  a  TA  phonon  replica  of  that  line.  The  problem 
is  that  the  1.513  line  is  present  in  the  unlmplanted  samples  and  the 
1.504  Line  is  not.  The  line  lias  not  been  reported  in  the  other  studies 
of  Mg-doped  or  implanted  CaAs.  When  the  J.49  eV  complex  was  resolved  on 
the  Dupont  310  lineshape  simulator,  the  L.504  line  was  found  to  be  fairly 
broad  compared  to  the  other  Jines.  There  is  no  guarantee  however  of  the 
uniqueness  of  the  fit  obtained  from  the  simulator  in  resolving  the  line- 
shape  into  its  component  lines. 

Tile  J  .5133  eV  line  has  been  identified  by  Yu  and  Park^^  and  Bogar- 
(92) 

dus  and  Bebb  as  an  impurity  exi  itou  complex.  in  particular,  it  is 

(92) 

associated  with  some  unknown  donor  or  donors  .  At  one  time  this  line 

was  thought  to  come  from  a  f ree  -exc i ton  emission.  The  energy  of  the  line 

(92) 

is  closest  to  the  exciton-ioni /ed  donor  line  in  Bogardus  and  Bebb 
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Beam  Current 


They  also  indicate  that  this  line  might  come  from  a  free  hole  to  neutral 
donor  recombination. 

These  are  only  suggested  assignments.  It  is  felt  that  the  free  to 
bound  line  for  Mg  at  L.492  eV  is  definite.  The  Zn  line  at  1.490  eV,  the 
Si  line  at  1.48  eV  and  their  phonon  replicas  are  also  reasonably  certain. 
The  others  are  less  so.  For  this  study  the  only  line  of  major  importance 
is  the  Mg  line.  While  the  copper  line  is  used  for  normalizing,  it  is  not 
important  to  know  what  caused  it  as  long  as  it  is  known  that  the  source 
of  this  line  is  uniformly  distributed  over  the  depth  of  the  epilayer. 

Spectral  Variation  with  Current  Density 

In  this  section  it  will  be  illustrated  how  the  spectrum  changes  as 
the  current  density  on  the  sample  varies.  The  current  density  is  speci¬ 
fied  in  terms  of  the  total  current  collected  by  the  Faraday  cup.  This  is 
the  quantity  that  is  measured.  As  the  cross  sectional  area  of  the  elec¬ 
tron  beam  is  held  constant  (by  virtue  of  the  cross  section  of  the  beam 
always  being  bigger  than  the  aperture  for  the  Faraday  cup),  constant 
total  current  implies  constant  current  densities,  assuming  the  beam  has  a 
uniform  crosssection.  Since  a  cross  section  of  the  beam  suffers  some 
nonuniformities,  the  following  procedure  was  adopted  to  minimize  its 
effect.  The  electron  beam  is  first  steered  onto  a  Faraday  cup  into  a 
position  that  maximizes  the  beam  current  into  the  cup.  Then  the  beam  is 
switched  to  t lie  sample  and  Che  amount  of  luminescence  at  a  particular 
line  is  maximized  by  moving  the  beam  around  on  the  sample.  This  proce¬ 
dure  tends  to  place  the  same  part  of  the  beam  on  the  sample  that  was 
measured  in  the  cup  and  to  locate  the  beam  on  the  sample  in  the  same 
place  as  was  used  for  the  previous  measurement.  Figure  (19)  is  the 
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luminescence  from  a  90  keV  Mg  implant  inLo  epitaxial  GaAs  at  a  fluence  of 
-2 

1EJ3  cm  .  The  spectra  were  taken  with  5.0,  1.0  and  0.1  microamp  beam 
currents  and  a  15  keV  electron  beam  energy. 

The  first  plot,  figure  (19-a),  is  of  the  spectrum  taken  with  a  5 
microamp  beam  current.  Figure  (19-b)  expands  the  spectrum  around  the 
1.492  and  1.455  eV  lines.  The  Mg  line  at  1.492  eV  is  the  major  line  in 
that  spectrum.  The  line  at  1.504  eV  is  missing  in  all  three  spectra. 

The  peak  at  1.455  eV  is  the  major  peak  in  the  phonon  replica  structure. 
This  supports  the  earlier  identification  of  the  1.455  eV  line  as  a  phonon 
replica  of  the  1.492  line.  As  the  beam  current  is  decreased  to  1.0 
microamp,  (figures  (19-c,  d)),  the  amplitude  jf  the  1.492  line  drops 
relative  to  the  other  lines.  In  fact,  the  Zn  line  is  almost  as  large  as 
the  Mg  line.  Simultaneously,  the  1.455  eV  line  drops  relative  to  the 
1.449  line.  A  new  line  at  1.453  eV  appears  that  is  37  meV  below  the  1.490 
eV  line.  This  line  is  a  phonon  replica  of  the  zinc  line.  It  could  not 
been  seen  before  because  the  zinc  line  was  relatively  too  small.  The 
third  spectrum  is  at  0.1  inicroainp.  Now  the  Mg  line  is  smaller  than  the 
Zn  line.  Similarly,  the  intensity  of  the  1.455  line  has  dropped  below 
the  1.453  eV  line.  The  llncsliapes  also  appear  to  be  more  clearly  resolved 
at  the  lower  currents. 

The  reason  for  the  dramatic  change  with  current  density  is  not  at 
ail  clear.  The  trends  were  the  same  in  ail  of  the  samples  examined.  A 
possible  explanation  is  that  the  probability  for  recombination  is  higher 
for  the  deeper  Lying  acceptors,  that  is,  the  probability  for  the  free  to 
bound  transition  involving  Zn  is  more  probable  tnan  the  same  transition 
with  Mg.  As  the  current  density  increases,  the  recombination  patli  through 
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WAVELENGTH 

Expanded  5  microamp 


RELATIVE  MAGNITUDE  (ARB  UNITS) 


Kif>.  19  (c)  !  mlcnwimp  lu-ain  current 
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since  there  are  less  Zn  impurities  Llian  Mg. 


Zn  saturates  belore  the  Mg, 

At  tills  point  ,  the  Mg  line  takes  over  ami  becomes  the  more  important 

(ij  5) 

line.  According  to  Ihimke's  theoretical  approach  for  the  recombina¬ 

tion  time,  t,  of  electrons  with  acceptor-bound  holes  at  Low  temperatures, 

t  u  L^/l>  (Z4) 

A  A 

E  is  the  acceptor  binding  energy  and  1’  the  concentration  of  bound 
A  A 

holes.  Since  E,  is  larger  for  Z n  than  Mg,  i’  must  also  be  larger  for  Zn 
A  A 

hy  enough  to  make  T  for  Zn  smaller  Llian  r  for  Mg.  What  happens  is  that 

the  holes  bound  to  the  impurities  seek  out  the  higher  impurity  line  in 

the  energy  gap  (Zn  over  Mg)  as  this  is  the  state  of  Lower  energy  for 

holes.  Thus ,  if  the  density  of  Zu  impurities  is  large  enough  compared  to 

the  Mg  density,  the  Zn  free  to  bound  transition  has  a  shorter  lifetime 

than  the  Mg  transition  due  to  the  larger  concentration  of  bound  holes. 

As  the  beam  current  increases  there  are  more  holes  available  to  be  bound. 

Eventually  the  Zn  recombination  path  is  saturated  as  the  number  of  bound 

holes  approaches  the  numbet  of  Zn  impurities.  As  there  are  mure  Mg  ions 

than  Zn  ions  available  to  bind  the  holes,  s'  for  Mg  starts  to  increase 

A 

and  Lin-  I  i  let  i me.  lor  Mg  decreases  to  the  point  where  Mg  is  the  more 
important  I  in-  .  A  similar  effect  should  hold  tor  the  Zn  and  Si  lines. 

As  the  current  density  goes  down,  tin  Si  1  ine  intensity  should  increase 
.umpired  to  t  he  Zn  line.  Table  V  shows  the  results  for  the  unimplanted/ 
i  i  i .  i  n  i  ■  •  -i  I  e;|  s  ample.  This  sample  was  -  hoseii  so  that  the  presence  of  tin-  Mg 
line  does  not  alft  t  the  results.  The  data  is  taken  at  lu  kcV .  l"ne  Si 


i  Dt*  i  i 

i  u<  teas  i  ng  t  as  t  e  i 

1  ban 

t  he  ..a  it  in-  as 

tin-  current  iner 

cases . 

I'll 

fi  n\y,i’ 

is  more  pronounced 

1  rom 

1.0  to  0.1  mil  t 

mumps  since  the 

larger 

Zn 
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peak  overshadows  Che  Si  peak  at  5  microamps.  The  ratio  is  based  on  the 
peak  of  the  line  as  measured.  The  peaks  do  not  have  the  effects  of  the 
other  lines  subtracted  out  before  the  ratio  is  calculated.  To  do  so 
would  make  the  changes  more  dramatic.  Other  samples  show  similar  trends. 
The  spectrum  shows  in  figure  (19)  shows  a  similar  trend  in  the  ratio  from 
1.0  to  0.1  microamp.  The  trend  is  reversed  from  1.0  to  5.0  microamps. 

This  reversal  is  caused  by  the  Si  and  Zn  lines  being  swamped  by  the  much 
larger  Mg  line.  The  reason  the  ratios  at  1.0  and  0.1  microamp  follow  the 
proper  trend  is  that  the  lines  are  better  resolved  at  these  beam  currents, 
hence  the  Mg  line  has  leas  effect. 

Table  V 

Ratio  of  Peak  Intensity  of  Si  to  Zn  Free 

Beam  Current  Unimplanted/ 

(microamp)  Unannealed  Sample 

5.0  .26 

1.0  .29 

0.1  .58 

*Spectrum  shown  in  figure  (19) . 

Another  point  to  be  made  is  that  the  peaks  are  better  resolved  at 
the  lower  current  densities.  This  is  attributed  to  less  local  heating  of 
the  sample  at  the  lower  electron  beam  current  densities.  Also  the  abso¬ 
lute  intensities  of  all  lines  decrease  as  the  current  density  drops 
because  there  are  less  electroncs  and  holes  available  for  recombination. 
Kurbatov  et  ajd04)  also  studied  the  influence  of  the  excitation 

level  on  GaAs  luminescence.  Their  data  was  taken  at  60  keV  beam  energy 

2 

with  current  densities  from  0.6  to  20  amps/cm  .  This  is  four  tiroes  the 


-to-Bound  Transition 

90  keV/lE13  cm"2 
Mg/GaAs  Sample* 

.32 

.29 

.  46 
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maximum  beam  energy  used  here  and  about  1000  to  320,000  times  the  current 
density.  Thus  their  results  are  not  directly  related  to  those  reported 
here.  Still,  they  observed  new  lines  emerging  as  the  beam  current  den¬ 
sity  increased.  These  lines  tended  to  be  associated  with  different  types 
of  transitions  instead  of  the  same  type  of  transition  through  different 
impurities.  For  instance,  at  10°K  in  n-type  GaAs  the  line  dominant  at 
low  current  densities  is  a  free-to-bound  transition  through  an  acceptor, 
while  the  dominant  line  at  higher  current  densities  appears  to  be  an 
exciton  bound  to  a  donor.  For  the  data  in  this  report,  the  type  of 
transition  is  a  free  electron  to  a  hole  bound  on  an  acceptor;  the  dif¬ 
ferent  lines  are  associated  with  different  acceptor  impurities. 

The  reason  that  the  variation  with  current  density  is  Important  to 
profiling  is  that  the  profile  information  is  contained  in  the  variation 
of  the  Mg  line  with  beam  energy.  The  higher  current  densities  are 
required  in  order  to  force  the  Mg  line  to  be  the  major  contributor  to  the 
1.49  eV  complex.  If  it  is  not,  the  profile  information  is  further  ob¬ 
scured  by  the  adjacent  lines.  This  results  in  fewer  details  of  the 
profiles  being  resolved  and  reduces  the  confidence  factor  in  the  final 
result.  While  a  larger  current  density  produces  a  larger  ratio  of  inten¬ 
sities  of  Mg  to  Zn  lines,  a  smaller  current  density  is  easier  to  achieve 
experimentally.  The  reason  for  this  is  that  at  the  lower  beam  energies 
the  electron  gun  becomes  space  charge  limited  and  the  maximum  gun  current 
is  reduced.  With  the  present  experimental  arrangement  5  microamps  at  2.5 
keV  is  the  best  that  could  be  done  and  only  occasionally  could  this  re¬ 
sult  be  reached.  Frequently  3  or  4  keV  at  5  microamps  is  the  best  that 
could  be  done.  The  L(V)  curves  were  determined  at  5  microamps  since  this 


seemed  the  best  compromise  between  low  energy  response  and  high  Intensity 
Mg  lines. 

Luminescence  versus  Beam  Voltage  Curves 

This  section  contains  the  primary  results  of  the  investigation.  The 
goal  of  this  research  is  to  determine  the  feasibility  of  using  cathodo- 
luminescence  as  a  depth  probing  tool.  In  this  section  the  experimental 
results  are  described  and  compared  with  the  theoretical  results  of  the 
last  chapter. 

First,  it  is  necessary  to  describe  again  how  the  data  was  taken  and 

presented.  As  noted  above,  the  beam  current  was  kept  at  5  microamps  when 

the  L(V)  curves  were  to  be  determined  using  the  Mg  line.  In  order  to 

cancel  any  variations  from  run  to  run,  such  as  changes  in  the  optics, 

changes  in  beam  position,  etc.,  the  Mg  line  was  normalized  to  the  copper 

line.  Since  the  calculation  predicts  the  total  number  of  photons  per 

second  to  be  observed,  it  is  necessary  to  determine  the  integrated  line 

intensity.  If  the  lineshapes  are  Gaussian, 

2 

-ax 

f(x)  -  c  J  6  (25) 

then  the  integrated  line  intensity  is 

2 

„  -ax 

A  -  /  c  -  e  dX  -  c  (26) 

-CO  ir 

(94) 

If  the  lineshapes  are  of  the  form  suggested  by  Eagles  for  band  to 
acceptor  recombination, 

f(x)  -  cx^exp  ax  (27) 

where  c  and  a  are  constants,  then  the  integrated  line  intensity  is, 
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(28) 


<o  -ax 

A  *  /  cx  e  dX 
o 

In  either  case,  the  peak  of  the  line  should  be  proportional  to  the  inte¬ 
grated  line  intensity  provided  that  the  lineshape  is  either  of  the  form 
of  equations  (25)  or  (27),  that  "a"  is  not  a  function  of  beam  energy,  and 
that  the  lines  are  clearly  resolved.  If  "a"  in  equations  (25)  and  (27) 
is  a  function  of  beam  energy  then  the  proper  way  to  scale  the  peak  line 
intensity  to  the  integrated  line  intensity  is 

AREA  a  PEAK/ /a  <29> 

for  a  Gaussian  lineshape  and 

AREA  a  PEAK/-  (30> 

a 

for  Eagle's  lineshape.  It  will  be  shortly  demonstrated  that  both  the 
peak  intensities  and  the  Integrated  lineshapes  yield  approximately  the 
same  L(V)  curves. 

L(V)  curves  determined  from  the  peak  and  area  of  the  lineshapes  are 
compared  in  figures  (20)  and  (21).  The  data  points  for  the  area  curves 
are  determined  by  resolving  the  Mg  line  from  the  1.49  eV  complex  using 
the  Dupont  310  lineshape  simulator.  The  Mg  integrated  line  intensity  can 
then  be  read  from  the  machine.  The  area  of  the  copper  line  can  also  be 
determined  using  the  simulator.  A  point  on  the  L(V)  curve  is  determined 
by  dividing  the  area  of  the  Mg  line  by  the  area  of  the  copper  line  and 
scaling  the  result  so  that  the  peak  of  the  L(V)  curve  is  at  400* 

The  peak  L(V)  curve  is  calculated  by  determining  "a"  from  the  data 
and  scaling  the  peak  according  to  equation  (29).  The  Gaussian  scaling 
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cr  2  ■ 886C 

3  "  3 

a-r  ott 


factor  was  chosen  since  Gaussian  curves  were  used  to  determine  the  area 

L(V)  curves.  A  data  point  on  the  L(V)  curve  is  calculated  by  dividing  the 

Mg  line  by  the  Cu  line  and  scaling  to  400  as  before.  The  results  for  the 

-2 

60  keV  implant  at  5E13  cm  dose  are  presented  in  Table  VI.  The  ratios 
have  been  scaled  to  400  for  the  largest  ratio  in  each  column.  Table  6  is 
plotted  in  figure  (20).  The  curves  are  very  close  considering  the  arbi¬ 
trariness  with  which  the  area  of  the  Mg  line  is  determined.  Figure  (21) 

-2 

shows  the  results  for  the  90  keV  implant  at  5E13  cm  .  While  the  results 
are  not  as  close  as  before,  the  same  trends  are  evident.  Thus,  although 
the  integrated  line  intensity  plot  is  necessary  for  quantitative  agree¬ 
ment,  the  peak  intensity  curves  are  at  least  qualitatively  right.  The 
curves  are  not  expected  to  agree  exactly  since  the  peak  intensity  curves 
are  determined  from  the  peaks  that  include  the  effects  of  all  the  competing 
lines  in  the  1.490  eV  complex. 


Table  VI 


Beam  Energy  (keV) 

Peak  Ratio 

Area  Ratio 

4 

400 

400 

5 

292 

295 

7 

267 

273 

9 

201 

212 

11 

181 

195 

13 

174 

177 

15 

130 

138 

Comparison  of  L(V)  Curves  for  Peak  vs  Integrated  Line- 
shape  Techniques  Mg  implanted  into  GaAs  at  60  keV/5E13  cm 


Electron  beam  energyckevj 


Fig.  20  Comparison  of  60  keV/5F13  Mg  Implanted  GaAs 


Mg  Implanted  GaAs 


Several  peak  Intensity  curves  were  obtained  In  order  to  get  a  quail' 
tatlve  feel  for  how  the  L(V)  curves  might  vary.  The  reason  for  obtaining 
peak  curves  instead  of  the  integrated  line  intensity  curves  is  that  they 
can  be  obtained  much  faster. 

Figures  (22)  and  (23)  are  peak  L(V)  curves  for  5E13  and  1E14  cm  ^ 

_2 

dose  Mg  implanted  GaAs.  The  peak  curves  for  5E13  cm  show  the  sort  of 

trends  expected.  The  120  keV  sample  yields  more  luminescence  from  higher 

energy  beams  than  does  the  60  or  90  keV  sample.  The  90  keV  curve  comes 

next  and  then  the  60  keV  curve.  While  the  curves  do  not  look  exactly 

like  the  curves  in  figure  (11-d)  the  same  general  trend  is  evident.  All 

the  curves  show  the  type  of  fall  off  predicted  in  figure  (11).  There  are 

_2 

two  apparent  problems  with  the  5E13  cm  curves.  The  first  one  is  that 
the  90  keV  curve  falls  off  much  faster  than  the  60  or  120  keV  curves. 

The  other  is  that  the  60  keV  curve  has  somewhat  of  a  peak  in  the  tail  of 

the  curve.  These  problems  will  be  addressed  in  the  next  section. 

-2  -< 
The  1E14  cm  curves  in  figure  (23)  are  very  similar  to  the  5E13  cm 

curves.  The  120  keV  curve  shows  more  luminescence  coming  from  the  higher 

energy  beams  thus  reflecting  the  fact  that  more  implanted  ions  exist  at 

deeper  depths  than  are  reached  with  the  90  keV  beam.  Again  the  same  type 

of  fall  off  is  seen,  although  the  120  keV  curve  appears  to  turn  up 

slightly  at  the  end. 

Earlier,  it  was  shown  that  the  higher  current  densities  enhance  the 
Mg  free-to-bound  line.  By  the  same  token,  lower  current  densities  favor 
the  Zn  1  Ine.  Figure  (24)  shows  a  comparison  of  the  L(V)  curves  deter¬ 
mined  from  the  peak  of  the  Mg  line  at  5  microamps  with  the  L(V)  curve 
from  the  peak  of  the  Zn  line  at  0.1  microamps.  The  Mg  line  clearly  shows 
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Fig.  24  Peak  L(V)  Curves  1  14  cm  /120  keV  Mg  Implanted  GaAs. 
Mg  Line  at  5  microamps  Beam  Current  and  Zn  Line  at 
0.1  microamps. 


the  effect  of  the  implant,  while  the  Zn  line  is  essentially  flat.  This 
behavior  would  suggest  that  the  Zn  is  uniform  with  depth,  as  the  ratio  of 
the  intensities  of  the  lines  from  two  uniformly  doped  impurities  should 
be  a  constant.  Since  Zn  was  present  in  the  epilayer  before  implantation, 
this  seems  like  a  valid  result. 

In  figure  (25)  the  computer  simulation  is  presented  for  two  uniform 

implants  with  a  diffusion  length  of  .21  microns,  a  surface  recombination 

velocity  of  1.65  x  10^  cm/ sec  and  absorption  coefficients  for  Zn  and  Cu 

of  70  and  1.5  cm  ^  respectively.  The  results  are  reasonably  close. 

While  the  data  points  bend  slightly  up  for  the  Zn  line,  the  simulation  is 

flat.  The  increase  may  reflect  a  slight  increase  in  the  Zn  density  with 

depth  (less  than  9%  change),  the  uncertainty  in  the  data  points,  or  a 

superl inear  variation  of  cathodoluminescence  intensity  with  net  carrier 
(191 

generation  rate  .  In  any  case,  there  is  a  clear  distinction  between 
the  appearance  of  the  Zn  line  and  the  Mg  line  both  in  theory  and  experiment. 

Figure  (26)  presents  the  three  L(V)  curves  considered  the  most 
likely  to  give  quantitative  agreement  with  theory.  These  were  the  last 
ones  taken  after  most  of  the  experimental  problems  had  been  worked  out. 

Many  of  the  same  observations  made  for  the  peak  L(V)  curves  apply  here 
since  the  peak  and  integrated  lineshape  curves  agree  quite  well.  The 
peak  in  the  tail  of  the  60  keV  curve  as  well  as  the  rapid  fall  off  of  the 
90  keV  curve  will  be  addressed  in  the  next  section  where  the  model  and 
the  data  will  be  put  together  to  see  how  well  they  agree. 

Reliability  of  Results 

Depth  profiles  were  made  by  Y.  K.  Yeo  of  the  U.S.A.F.  Avionics 

_o 

Laboratory  on  the  5E13  cm  /90  and  120  keV  implants.  The  depth  profiles 
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Fig.  25  Comparison  of  Calculated  and  Experimental  L(V)  Curves 
for  Two  Uniform  Impurities.  Dots  are  experimental 
points . 
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Fig.  26  Integrated  Lineshape  L(V)  Curves  for 
60,  90,  120  keV  Mg  Implanted  GaAs  at 
5E1  3  i:m“^  Fluence 


showed  a  considerable  scatter  in  the  data  points,  although  they  did 
vaguely  resemble  profiles  determined  by  him  earlier  for  other  Mg  im¬ 
planted  samples.  There  was  a  discoloration  left  on  the  samples  after 
the  L(V)  curves  had  been  obtained.  In  fact,  as  the  GaAs  was  etched  away 
from  one  sample,  the  discoloration  actually  seemed  to  get  bigger.  It  had 
been  supposed  that  this  was  a  thin  film  of  carbon  resulting  from  the 
interaction  of  the  electron  beam  and  the  hydrocarbons  from  possible 
backdif fusion  of  the  diffusion  pump  oil  or  the  aquadag  that  covered  the 
copper  shield  that  holds  the  samples  in  place.  Further  inquiry  then 
turned  up  the  fact  that  the  etch  would  not  dissolve  carbon.  A  surface 
height  measurement  was  made  with  a  Dektak  stylus  which  confirmed  that  a  2 
to  5  micron  layer  had  built  up  on  the  surface.  As  the  sample  was  etched, 
the  GaAs  not  protected  by  the  film  was  etched  away.  This  leads  to  the 
appearance  of  profiling,  but  does  not  accurately  reflect  the  impurity 
profile.  This  explains  the  scatter  in  the  data  points  as  different  parts 
of  the  surface  are  etched  at  different  times  and  rates.  It  also  explains 
the  apparent  increase  in  the  spot  size  because  of  the  shadowing  effect  of 
the  film  protective  layer  as  the  GaAs  is  etched  away  around  it. 

This  brings  up  the  question  of  how  the  L(V)  curves  are  affected  by 
the  film  layer.  First,  a  carbon  layer  does  not  noticeably  attenuate  the 
incident  electrons^^ .  What  it  does  to  the  emerging  luminescence  is 
another  question.  Fortunately,  since  the  L(V)  curve  is  determined  by 
dividing  the  Mg  by  the  Cu  line,  the  only  assumption  needed  to  ignore  the 
absorption  effect  of  the  carbon  is  that  the  absorption  coefficients  for 
luminescence  from  the  Mg  and  Cu  impurities  have  the  same  functional 
dependence  with  film  thickness.  If  they  do,  then  taking  the  ratio  of  two 


lines  cancels  out  any  effects  of  absorption  from  the  film.  Nagy  et 

^2.(96)  show  the  transmittance  of  a  carbon  film  in  the  8000  to  10,000  A 

range  as  smoothly  varying  from  50  to  60%.  It  seems  unlikely  that  the 

variation  of  the  absorption  coefficients  with  film  thickness  will  be 

significantly  different  for  Mg  and  Cu.  Therefore,  the  L(V)  curves  are 

assumed  unaffected  by  the  possible  existence  of  a  carbon  film. 

Figure  (27)  shows  the  calculated  and  experimental  L(V)  curves  for 
_2 

the  120  keV/ 5E13  cm  sample.  All  things  considered  the  agreement  is 

remarkable.  The  values  used  in  the  calculation  are: 

L  =0.21  microns 
n 

ALPHA  w  =  70  cm  1 
Mg 

ALPHA  _  *1.5  cm-1 

Cu 

s  =  1.65  x  10°  cm/sec 

where  L  is  the  ambipolar  diffusion  length,  ALPHA  is  the  absorption  coef 
flcient  and  s  is  the  surface  recombination  velocity. 

The  question  arises  as  to  whether  these  are  reasonable  values  to 
use.  Values  at  10°K  are  hard  to  find.  The  diffusion  length  can  be  cal¬ 
culated  from 


L  - 


(31) 


where  k  is  Boltzmann’s  constant,  T  is  temperature,  e  the  electronic 
charge,  and  b  the  mobility.  Bolger^*^  et  al  measured  electron  mobi- 

4 

lities  down  to  4.2°K.  At  108K  they  found  mobilities  ranging  from  1  x  10 

to  4  x  10^  cm^/v  sec  for  donor  impurity  concentrations  from  4*9  x  10^  to 

3.0  x  10  cm  .  Radiative  recombination  times  for  Impurities  to  the 

valence  or  conduction  band  can  be  estimated  from  Dumke'e  theoretical 

17  -3 

result.  Assuming  the  number  of  acceptors  to  be  3  x  10  cm  which  is 
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consistent  with  Yeo's  measured  result  and  a  5L13  cm  dose  spread  over 

-9 

1.67  microns,  the  radiative  lifetime  at  10°K  is  6.9  x  10  seconds. 

Since  nonradiative  and  other  radiative  channels  shorten  this  figure  it 
has  to  be  considered  an  upper  limit.  Using  both  upper  limits  for  T  and 
mobility  and  calculating  the  diffusion  length  from  equation  (31)  leads  to 
the  conclusion  that  the  maximum  diffusion  length  should  be  0.28  microns. 
This  is  not  far  from  the  0.21  microns  used  in  the  computer  fit.  Since 
0.28  microns  is  considered  an  upper  limit  the  0.21  micron  figures  seems 
very  reasonable. 

The  absorption  coefficients  in  high  resistivity  GaAs  were  measured 

by  Sturge^^^  from  10°  to  294°K.  From  0.6  to  1.4  eV  the  absorption 

coefficient  is  small,  staying  under  4  cm  * .  As  the  band  gap  of  GaAs  is 

approached  the  absorption  coefficient  rises  rapidly  to  a  value  around  1  x 
4  -1 

10  cm  .  The  absorption  coefficient  at  10°K  around  1.49  eV  is  between  60 
and  80  cm  ^ .  These  values  are  consistent  with  the  values  measured  by 
Redfield  and  Afromowitz^®^  for  p  type  GaAs  at  5°K.  Thus  the  values 
used  in  the  computer  fit  for  the  absorption  coefficients  are  supported  by 
other  reported  values.  In  addition  the  results  are  not  particularly 
sensitive  to  changes  in  the  copper  absorption  coefficient.  It  can  shift 
by  an  order  of  magnitude  either  way  and  not  significantly  change  the 
results. 

The  surface  recombination  velocity  falls  right  on  the  curve  given  by 

Jastrzebskl  et  al^^^  for  similarly  doped  samples.  Their  value  for  a 

1 7  •“  3  6 

mid  10  cm  doped  sample  is  about  1.8  x  10  cm/sec.  This  result  is  for 

room  temperature  but  Vilms  and  Spicer  observed  that  the  surface  recombina¬ 
tion  velocity  did  not  change  much  over  the  80°  to  300°K  range.  Assuming 


that  the  surface  recombination  velocity  does  not  change  with  temperature 
then  this  value,  too,  is  considered  about  right. 

The  final  conclusion  concerning  the  values  chosen  for  the  computer 

fit  is  that  they  are  reasonable  based  on  other  workers  results  and  theo- 

-2 

retical  calculations.  The  two  other  5E13  cm  dose  curves  show  similar 
results. 

The  90  keV  curve  is  fitted  using  the  same  diffusion  length,  coeffi¬ 
cients  and  surface  recombination  velocity  as  the  120  keV  sample.  The 
results  are  shown  in  figure  (28).  The  fit  is  not  nearly  as  good  as 
before.  The  fit  can  be  improved  by  going  to  a  much  shorter  diffusion 
length  and  higher  surface  recombination  velocity.  A  shorter  diffusion 
length  can  be  supported  in  the  literature;  a  higher  surface  recombination 
velocity  Is,  however,  much  harder  to  justify.  A  good  computer  fit  to  the 
experimental  data  could  not  be  obtained  since  below  a  diffusion  length  of 
about  0.1  micron  the  simulation  failed  to  converge. 

The  60  keV  (.(V)  curve  is  fitted  assuming  a  slightly  different  impur¬ 
ity  profile.  This  profile  consists  of  two  LSS  profiles  added  together. 

The  first  is  the  LSS  profile  predicted  by  LSS  theory  for  a  60  keV  implant. 
The  second  profile  has  the  same  shape  as  the  first  but  is  centered  around 
a  0.5  micron  depth  Lnstead  of  the  0.0bl5  micron  depth  predicted  by  LSS 
and  has  an  amplitude  equal  to  0.55  times  the  first  peak.  With  this 
change  and  assuming  the  same  parameters  as  the  120  keV  implant  gives  the 
fit  shown  in  figure  (29).  Figure  ( 30)  shows  the  profile  used  in  figure 
(29).  The  physical  nature  ot  the  second  hump  is  unknown  but  could  be 
caused  by  a  component  of  the  implant  beam  penetrating  down  a  crystalline 
axis  (channeling) . 


The  experience  with  the  60  keV  implant  brings  up  a  question  as  to 
the  resolution  possibilities  of  this  technique  as  regards  double  peaks. 

In  order  to  address  this  question  a  study  was  undertaken  to  determine  the 
range  over  which  two  peaks  make  measurable  changes  in  the  L(V)  curves  for 
a  double  profile  as  is  found  in  the  60  keV  sample.  The  details  are  con¬ 
sidered  in  Appendix  B.  The  results  indicate  that  with  the  parameters 
used  here,  peaks  separated  from  0.25  to  slightly  less  than  10.  microns 
could  be  distinguished. 

An  important  point  to  note  in  figures  (27),  (28),  and  (29)  is  that 
data  points  below  4  keV  play  an  important  part  in  distinguishing  one  pro¬ 
file  from  another.  It  <s  anticipated  that  a  experimental  setup  capable 
of  yielding  data  points  down  to  election  beam  energies  of  2  keV  is  very 
important  in  identifying  different  implant  profiles. 

Summary  of  Important  Results 

This  section  is  a  summary  of  the  most  important  results  of  this 
study.  The  first  result  obtained  is  in  the  depth-dose  curves  for  elec¬ 
tron  beam  penetration  into  GaAs.  These  curves  are  important  for  both 
cathodolumlnescent  and  microprobe  work  as  they  help  to  quantify  models  of 
both  processes. 

The  next  result  Is  the  spectra  from  various  samples  of  Mg  Implanted 
epitaxial  GaAs.  These  spectra  had  not  been  studied  before.  Several 
lines  of  the  spectra  between  1.3  and  1.5  uV  were  identified.  The  effect 
on  the  spectra  of  changing  the  current  density  was  studied.  In  this 
section,  it  was  shown  how  dramatically  the  intensity  of  various  lines 
change  with  current  density.  A  proposed  explanation  of  this  variation 
was  given  for  the  free-to-bound  transitions. 
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29  Comparison  of  Calculated  and  Experimental 

I.(V)  Curves  for  60  keV/5E13  cm-2  Mg  Implanted 
OaAs.  Dots  are  experimental  points. 


line  for  a  Gaussian  lineshape  and  the  llneshape  proposed  by  Eagles. 

Next,  the  theoretical  and  experimental  curves  were  compared.  Good  agree¬ 
ment  is  found  for  these  curves  providing  the  proper  parameters  are  chosen. 
The  reliability  of  these  choices  is  discussed  and  it  is  concluded  that 
the  parameters  chosen  are  reasonable. 


V.  Conclusions  and  Recommendations 


The  purpose  of  this  research  is  to  "establish  a  quantitative  basis 
for  cathodoluminescence"  in  GaAs.  To  that  end  a  computer  model  was 
developed  of  the  process  and  data  was  taken  in  an  effort  to  validate  the 
model.  As  demonstrated  in  the  last  section  of  the  previous  chapter  the 
agreement  is  very  good.  The  conclusion  is  that  a  basis  has  been  estab¬ 
lished  for  quantitative  cathodoluminescence. 

There  are  several  other  conclusions  that  can  be  drawn  from  the 
results  of  this  study.  The  first  is  that  depth-resolved  cathodolumine¬ 
scence  can  be  used  not  only  to  identify  what  impurities  are  present  in 
GaAs,  but  also  where  the  concentration  is  highest.  This  was  demonstrated 
for  the  Si  impurity  where  it  was  shown  that  Si  diffused  into  the  epitax¬ 
ial  layer  from  the  surface  and  not  from  the  substrate. 

Next,  a  relation  between  the  diffusion  length,  recombination  time, 
diffusion  constant  and  surface  recombination  velocity  was  found.  This 
results  in  a  great  reduction  in  the  number  of  combinations  of  parameters 
that  need  be  examined  in  order  to  fit  experimental  and  theoretical  curves. 
It  also  means  that  without  some  independent  determination  of  one  of  the 
parameters,  they  cannot  be  absolutely  determined  from  cathodoluminescence 
data. 

The  study  of  the  spectra  from  several  samples  failed  to  turn  up  any 
donor-acceptor  pair  transition.  This  supports  the  contention  of  Bishop 
et  al^^)  wh0  interpreted  Yu  and  Park's  data  as  an  implant  depth  gradient 
effect  and  not  a  donor-acceptor  pair  transition. 

In  the  section  of  the  effect  of  changing  the  current  density  on  the 
spectra,  the  dramatic  changes  in  the  spectra  that  can  take  place  were 
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shown.  Unless  some  care  is  taken  to  make  sure  that  the  current  density 
is  held  constant,  conclusions  drawn  from  spectral  changes  must  be  ques¬ 
tioned.  Setting  the  total  beam  current  to  a  constant  value  is  not  as 
important  as  maintaining  a  constant  current  density  at  the  surface. 

By  implication,  the  diffusion  lengths  for  annealed  Mg  implanted  GaAs 
at  10°K  must  be  very  short.  The  best  fit  for  the  luminescence  data  was 
at  a  .21  micron  diffusion  length,  much  shorter  than  that  suggested  in 
reference  106.  This  short  diffusion  length  means  the  excited  carriers 
recombine  shortly  after  excitation.  The  volume  of  GaAs  in  which  radia¬ 
tive  recombination  occurs  is  thus  not  much  greater  than  the  volume  ex¬ 
cited  by  the  electron  beam.  Thus,  the  spatial  resolution  of  cathodo- 
luminescence  is  about  that  of  the  electron  energy  loss  curves.  Profiling 
should  thus  be  possible  for  any  impurity  imp lent  profile  whose  dimensions 
are  on  the  order  of,  or  more  than  those  spanned  by  the  energy  loss  curves. 

The  recommendat ions  generally  concern  improvements  in  the  experimen¬ 
tal  facilities  and  theoretical  model.  They  are: 

(1)  improve  the  electron  gun  design  so  that  a  uniform  beam  can  be 
obtained  from  1  to  15  keV. 

(2)  Improve  the  gun  design  so  that  more  than  5  microamps  of  beam 
current  can  be  provided  at  1  keV  when  the  beam  is  defocused  to  cover  the 
entire  entrance  aperture  to  the  Faraday  cup. 

(3)  Include  the  effects  of  the  depth  dependence  of  the  diffusion 
length  In  the  calculation. 

(4)  Work  on  solving  the  Inverse  problem  of  deducing  the  profile, 
given  the  L(V)  curve. 

(5)  Study  more  samples  under  different  conditions  to  establish  a 
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larger  data  base  on  which  to  draw  conclusions. 

Of  these  recommendations,  numbers  (4)  and  (5)  are  the  most  important. 
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Appendix  B 


In  order  to  fit  the  computer  model  profile  to  the  60  keV/5E13  cm 
Mg  implanted  GaAs  sample  it  was  necessary  to  add  two  LSS  profiles  to¬ 
gether.  This  brings  up  the  question  as  to  how  well  can  two  peaks  be 
resolved  by  measuring  the  L(V)  curve.  In  order  to  give  some  insight  into 
this  problem  a  series  of  L(V)  curves  for  the  60  keV  implant  were  calcu¬ 
lated.  For  these  curves  the  same  values  of  diffusion  length,  surface 
recombination  velocity  and  absorption  coefficients  are  used  as  previously 
used  for  the  60  keV  Implant. 

Figure  (B-l)  shows  the  results  as  the  first  Caussian  is  held  at 
0.0615  micron  and  the  second  Gaussian  is  moved  from  0.0615  micron  to  3.0 

microns.  The  Gaussians  both  have  the  same  amplitude  and  standard  devia¬ 

tion  as  the  LSS,  60  keV  Mg  implant.  Two  curves  are  shown  in  each  plot. 
Where  the  two  can  be  distinguished  the  upper  curve  is  the  result  of  two 
Gaussians.  The  lower  curve  is  the  L(V)  curve  resulting  from  a  60  keV, 

LSS  implant.  Thus  the  effect  of  the  second  peak  can  be  compared  to  the 
LSS  profile  in  each  case. 

In  the  initial  plot  the  two  Gaussian  implant  profiles  lie  on  top  of 
one  another.  Since  there  is  an  arbitrary  amplitude  scaling  done  on  the 
L(V)  curve  the  two  L(V)  curves  are  identical.  As  the  second  peak  moves 
farther  into  the  crystal  the  L(V)  curve  begins  to  show  the  effect  of  the 

second  peak.  With  as  little  as  .0635  micron  separation  of  the  peaks  a 

significant  change  in  the  calculated  L(V)  curve  results.  There  seems  to 
be  an  optimal  depth  at  which  the  two  peaks  can  be  resolved.  For  the 
values  used  here  it  is  about  0.5  of  a  micron  depth  for  the  second  peak. 

As  the  second  peak  penetrates  even  deeper  into  the  GaAs,  the  impact  of 
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B-l  L(V)  Curves  for  Twin  Peak  Profiles.  RP1  -  Depth  of  Second  Peak 
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the  second  peak  diminishes.  This  is  because  less  of  the  electron  beam 
penetrates  to  this  depth  and  more  of  the  luminescence  is  absorbed  before 
it  reaches  the  surface.  In  the  limit  of  3.0  microns  there  is  no  effect 
of  the  second  peak  reflected  in  the  L(V)  curve. 

While  only  one  specific  set  of  parameters  have  been  considered  here 
two  general  conclusions  are  evident.  First,  there  is  little  or  no 
effect  on  the  L(V)  curve  for  very  closely  spaced  or  very  widely  spaced 
peaks.  Second,  there  is  an  optimal  separation  of  the  first  and  second 
peaks  at  which  the  L(V)  shows  tiie  biggest  change  compared  to  the  single 
LSS  profile. 
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